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O ■ Abstract 

■^ ' Let G = (y, E) be a graph with / : V — >■ Z_|_ a function assigning degree bounds to vertices. 

tH- . We present the first efficient algebraic algorithm to find an /-factor. The time is 0(/(F)"). 

CS| ' More generally for graphs with integral edge weights of maximum absolute value W we find 

a maximum weight /-factor in time 0{Wf{V)'^). (The algorithms are randomized, correct 

r/^' with high probability and Las Vegas; the time bound is worst-case.) We also present three 

specializations of these algorithms: For maximum weight perfect /-matching the algorithm is 
considerably simpler (and almost identical to its special case of ordinary weighted matching). 
For the single-source shortest-path problem in undirected graphs with conservative edge weights, 
we present a generalization of the shortest-path tree, and we compute it in 0{Wn'^) time. For 
bipartite graphs, we improve the known complexity bounds for vertex capacitated max-flow and 
min-cost max-flow on a subclass of graphs. 
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(^ ' 1 Introduction 

\q' 

Tij- I 6-matching and /-factors are basic combinatorial notions that generalize non-bipartite matching, 

CD ■ min-cost network flow, and others. This paper presents the first efficient algebraic algorithms for 

both weighted and unweighted 6-matchings and /-factors. Our algorithms for this broad class of 

problems are the most efficient algorithms known for a subclass of instances (graphs of high density, 

low degree constraints and low edge weights). We also discuss single-source all-sinks shortest paths 

r> ' in conservative undirected graphs. (There is no known reduction to directed graphs.) We prove 

C^ , the existence of a simple shortest-path "tree" for this setting. We also give efficient algorithms - 

combinatoric for sparse graphs and algebraic for dense - to construct it. 

We must ffist define 6-matching and /-factors. The literature is inconsistent but in essence we 
follow the classification of Schrijver |27) . For an undirected multigraph G = {V,E) with a function 
/ : V — 7- Z4., an f -factor is a subset of edges wherein each vertex v £ V has degree exactly f{v). 
For an undirected graph G = (y,E) with a function b : V ^ Z-|_, a (perfect) b-matching is a 
function x : E ^ Z+ such that each v GV has YlwvweE •^i'^''^) ~ ^(^)- ^^^^ ^^^^ ^'^^^ 6-matchings 
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have an unlimited number of copies of each edge makes them decidedly simpler. For instance b- 
matchings have essentially the same blossom structure (and linear programming dual variables) as 
ordinary matching |27| Ch.31]. Similarly our algorithm for weighted 6-matching is almost identical 
to its specialization to ordinary matching (6 = 1). In contrast the blossoms and dual variables for 
weighted /-factors are more involved |271 Ch.32] and our algorithm is more intricate. Thus our 
terminology reflects the difference in complexity of the two notions|l| 

The paper begins with unweighted /-factors, i.e., we wish to find an /-factor or show none exists. 
Let (/) = f(y) (or b{V)). We extend the Tutte matrix from matching to /-factors, i.e., we present a 
(f) X (f) matrix that is symbolically nonsingular iff the graph has an /-factor. Such a matrix can be 
derived by applying the Tutte matrix to an enlarged version of the given graph, or by specializing 
Lovasz's matrix for matroid parity |18| . But neither approach is compact enough to achieve our 
time bounds G Then we reuse the elimination framework for maximum cardinality matching, due to 
Mucha and Sankowski |21| and Harvey |17| . This allows us to find an /-factor in 0{(f)'^) randomized 
timelfl For dense graphs and small degree-constraints this improves the best-known time bound of 
0{\/^m) |12| . although the latter is deterministic. 



Time 


Author 


Oin'B) 


PuUeyblank (1973) [23] 


0{n^mlogB) 


Marsh (1979) [20] 


0{m'^ log n log B) 


Gabow (1983) [H 


[n'^'m + n log B{m + n log n)) 


Anstee (1987) [2] 


0{n^ log n{m + n log n)) 


Anstee (1987) [2] 


0{W(l)^) 


this paper 



Table 1: Time bounds for maximum 6-matching. B denotes max„6(f). 



Complexity 


Author 


0{(l)n-^) 


Urquhart (1965) [32] 


0{(j){in + nlogn)) 


Gabow (1983+1990) [121 HI] 


0(W(l)'') 


this paper 



Table 2: Time bounds for maximum weight /-factors on simple graphs. 



We turn to the more difficult weighted version of the problem. Here every edge has a numeric 
weight; for complexity results we assume weights are integers of magnitude < W . We seek a 
maximum f -factor, i.e., an /-factor with the greatest possible total weight. Efficient algebraic 
algorithms have been given for maximum matching (/ = 1) in time 0{Wn^), ffist for bipartite 
graphs [26] and recently for general graphs 



^Another version of 6-matching considers h{v) as an upper bound on the desired degree of v. This easily reduces 
to weighted perfect 6-matching by taking 2 copies of G joined by zero-weight edges. On the other hand a capacitated 
h-matching is defined by giving an upper bound u(e) to each value x{e). The simplicity of the uncapacitated case is 
lost, and we are back to /-factors. 

^The Tutte matrix becomes too large, mx m. Lovasz's matrix is (f> x (p but can involve integers that are too large, 
size n" or more. Our matrix only involves integers ±1. 

'^0{n") is the time needed for a straight-line program to multiply two n x n matrices; the best-known bound on 
a; is < 2.3727 [33]. 

*The O notation ignores factors of log {n(f>W). 



The usual approach to generahzed matching problems is by problem reduction. For instance 
in |27| . Ch.31 proves the properties of the 6- matching linear program and polytope by reducing 
to ordinary matching via vertex splitting; then Ch.32 reduces /-factors (capacitated 6-matching) 
to 6-matching. Efficient algorithms also use vertex splitting |12| or reduction to the bipartite case 
(plus by further processing) [2]. But reductions may obscure some structure. To avoid this our 
algorithms use a direct approach, and we get the following rewards. For 6-matching, as mentioned, 
the similarity of blossoms to ordinary matching blossoms leads to an algorithm that is no more 
involved than ordinary matching. For undirected shortest paths we get a simple definition of a 
generalized shortest-path tree. (Again such a definition may have been overlooked due to reliance 
on reductions, see below.) For /-factors we get a detailed understanding of the more complicated 
versions of the structures that first emerge in 6-matchings (2-edge connected components giving the 
cyclic part of blossoms - see Section 18. ip and in shortest paths (bridges giving the incident edges 
of blossoms - these correspond to the (ungeneralized) shortest-path tree - see Section [8. 2p . 

All three of our non-bipartite algorithms are implementations of the "shrinking procedure" given 
in |14| (a variant is the basis of the weighted matching algorithm of [7]). This procedure gives a direct 
way to find the optimum blossoms for a weighted /-factor - simply put, each blossom is (a subgraph 
of) a maximum weight "2/-unifactor" (a type of 2/-factor) in the graph with (the cyclic part of) 
all heavier blossoms contracted (see |14| or Section [7]). Note that the classic weighted matching 
algorithm of Edmonds [8] finds the optimum blossoms, but only after forming and discarding various 
other blossoms. So this approach does not provide a direct definition of the optimum blossoms. 

The first step of our algorithms use our generalized Tutte matrix to find the optimum dual 
variables of the vertices. Then we execute the shrinking procedure to get the blossoms, their duals, 
and a "weighted blossom tree" that gives the structure of the optimum /-factor. (This step is 
combinatoric. It is based on the detailed structure of 2/-unifactors that we derive.) The last step 
finds the desired /-factor using a top-down traversal of the weighted blossom tree: At each node we 
find an /-factor of a corresponding graph, using our algorithm for unweighted /-factors. In summary 
our algorithms (like [7] for ordinary matching) can be viewed as a (combinatoric) reduction of the 
weighted /-factor problem into two subproblems: finding the optimum dual variables of the vertices, 
and finding an unweighted /-factor. 

To facilitate understanding of the general /-factor algorithm we begin by presenting its special- 
ization to two subcases. First 6-matching. The blossoms, and hence the 26-unifactors, differ little 
from ordinary graph matching. As a result our development for weighted 6-matching is essentially 
identical the special case of ordinary matching, in terms of both the underlying combinatorics and 
the algorithmic details. When specialized to ordinary matching our algorithm provides a simple 
alternative to [7] . In fact an advantage is that our algorithm is Las Vegas - the dual variables allow 
us to check if the 6-matching is truly optimum. (Our approach to weighted matching/6-matching 
differs from [^ - at the highest level, we work with critical graphs while [7] works with perfect 
graphs.) 

Next we discuss shortest paths in undirected graphs with a conservative weight function - 
negative edges are allowed but not negative cycles. The obvious reduction to a directed graph 
(replace undirected edge uv by directed edges uv, vu) introduces negative cycles, and it is unclear 
how to handle this problem by the usual shortest-path techniques. 

We consider the single-source all-sinks version of the problem. Again, this problem is often solved 
by reduction, first to the single-source single-sink version and then to perfect matching, using either 
T-joins [27^ pp. 485-486] or vertex-splitting [27^ p. 487]. A path can be viewed as a type of 2- factor. 



(For instance an a6-path is an /-factor if we enlarge G with a loop at every vertex v £ V and set 
f{v) = 2 for V £ V — {a, b}, f{a) = f{b) = 1.) This enables us to solve the all-sinks version directly. 
Examining the blossom structure enables us to define a generalized shortest-path tree that, like 
the standard shortest-path tree for directed graphs, specifies a shortest path to every vertex from 
a chosen source. It is a combination of the standard shortest-path tree and the blossom tree. We 
give a complete derivation of the existence of this shortest-path structure, as well as an algebraic 
algorithm to construct it in time 0{Wn'^). We also construct the structure with combinatoric 
algorithms, in time 0{n{m + nlogn)) or 0{^Jna{rn, n) log n mlog{nW)). These bounds are all 
within logarithmic factors of the best-known bounds for constructing the directed shortest-path tree 

[23 Ch. 8], [MIES]. 

Although the shortest-path problem is classic, our definition of this structure appears to be 
new. Most notably, Sebo has characterized the structure of single-source shortest paths in undirected 
graphs, first for graphs with ±1 edge weights |29) and then extending to general weights by reduction 
|30) . Equation (4.2) of |29| (for ibl-weights, plus its version achieved by reduction for arbitrary 
weights) characterizes the shortest paths from a fixed source in terms of how they enter and leave 
"level sets" determined by the distance function. [29| also shows that the distances from the source 
can be computed using 0{n) perfect matching computations. Our structure differs from |29| I30| : 
it does not give a necessary and sufficient condition to be a shortest path, but it gives an exact 
specification of a specific set of shortest paths that are simply related to one another (as in the 
standard shortest-path tree). Note that one can give an alternative proof of the existence of our 
structure by starting from the results of |29| I30| . 

The general algorithm for maximum /-factors is the most difficult part of the paper. It involves 
a detailed study of the properties of blossoms. A simple example of how these blossoms differ from 
ordinary matching is that the hallmark of Edmonds' blossom algorithm - "blossoms shrink" - is not 
quite true. In other words for ordinary matching a blossom can be contracted and it becomes just 
an ordinary vertex. For /-factors we can contract the "cycle" part of the blossom, but its incident 
edges remain in the graph and must be treated differently from ordinary edges (see Section 18. 2p . 
Our discussion of shortest paths introduces this difficulty in the simplest case - here a blossom has 
exactly 1 incident edge (as opposed to an arbitrary number). Even ignoring this issue, another 
difficulty is that there are three types of edges that behave differently (see |14| . or Lemmas [33] 
and [35|) and the type of an edge is unknown to the algorithm! Again the three types are seen 
to arise naturally in shortest paths. Our contribution is to develop the combinatoric properties 
of these edges and blossoms so the shrinking procedure can be executed efficiently, given only the 
information provided by the Tutte matrix in our algebraic algorithm. 

While non-bipartite graphs present the greatest technical challenge, we also achieve some best- 
known time bounds for two bipartite problems, maximum network flow and min-cost network flow. 
Bipartite /-factors generalize network flow: max-flow (min-cost max-flow) is a special case of un- 
weighted (weighted) bipartite /-factors, respectively e.g. |15| . The question of an efficient algebraic 
max-fiow algorithm has confronted the community for some time. The only advance is the algorithm 
of Cheung et. al. [5], which checks whether d units of flow can be sent across a unit-capacity network 
in 0{d^~^m) time. We consider networks with integral vertex and edge capacities bounded by D. 
We find a max-flow in time 0{{Dn)'^) time and a min-cost max-flow in 0{W{Dn)^) time. The 
latter algorithm handles convex edge cost functions (with integral break-points) in the same time 
bound. The max-flow problem has a rich history (see e.g. |27^ Chs. 10, 12]) and our time bounds 
are the best-known for dense graphs with moderately high vertex capacities. Specifically, previous 
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Table 3: Sections for each step of the algorithm. X is an algebraic algorithm, C is combinatoric. 
Shortest-path distances can be found algebraically or combinatorially. 



algorithms for vertex-capacitated max-flow in dense networks (i.e., m = 0(7i^)) use 0(?i^/logn) 
time [1] or 0{n^'^ log D) time |16| . Previous algorithms for dense graph min-cost max-flow use 
0(n^ log D) time |10| or 0(n^ log n) time |22| . Previous algorithms for minimum convex-cost max- 
flow use 0{Dn'^ log D) time (by simple reduction to min-cost max-flow) or in O {ri'^ log D log (nW)) 
time [T5] . 

In summary the novel aspects of this paper are: 

• new time bounds for the fundamental problems of 6-matching, undirected single-source shortest 
paths, and /-factors; 

• extension of the Tutte matrix for matching to /-factors; 

• definition of the shortest-path structure for undirected graphs, plus algebraic and combinatorial 
algorithms to construct it; 

• an algebraic algorithm for 6-matching that is no more involved than ordinary matching; 

• an algebraic algorithm for /-factors based on new combinatorial properties of blossoms, which 
are qualitatively different from matching blossoms; 

• new time bounds for vertex-capacitated max-flow, min-cost max-flow and minimum convex- 
cost max-flow on dense graphs. 



Organization of the paper The next two sections define our terminology and review the alge- 
braic tools that we use. 

Section |4] gives an overview of the entire paper, by discussing the special case of bipartite /- 
factors. (The reader should bear in mind that the non-bipartite case is the highlight of this paper. It 
is technically much more demanding.) In detail, Section |4] starts with an 0{(j)'^) time algorithm for 
unweighted bipartite /-factors. This requires our generalized Tutte matrix (Section 14. Ij Section [TT] 
extends this to non-bipartite graphs) plus the algorithmic details (Section l4.2| Section [T2l generalizes 
these details to non-bipartite graphs). Section 14.31 gives the algorithm for weighted bipartite /- 
factors. Flows are discussed in Section [4.41 



As mentioned above the algorithms for general graphs have three steps: 

(i) find the weights of factors of perturbed graphs; 

(a) use the weights to construct a blossom forest; 

(in) use an unweighted algorithm to extract an optimal solution from the blossom forest. 

These steps are implemented in different sections, depending on the problem of interest, as indicated 
in Table 3. For instance the complete 6-matching algorithm starts with the algorithm of Section [T3l 
to get the dual variables; then it constructs a weighted blossom forest and traverses the forest to 
extract the optimal solution, using algorithms in Section [5j the traversal uses the algorithm of 
Section [12] to find various unweighted 6-matchings. 

The most involved part of the paper is the construction of blossom forest for /-factors in Sec- 
tion [HI This section is preceded by two sections that introduce our combinatoric ideas in simpler 
settings: Section [5] gives the combinatoric portion of our weighted 6-matching algorithm. Section [6] 
proves the existence of our shortest-path structure. Section [7] reviews fundamental background ma- 
terial on /-factors (illustrating it by shortest-paths). After the combinatoric algorithm for weighted 
/-factors in Section [HI Section [9] gives the algorithm to construct the shortest-path structure, and 
Section [TOJ gives the combinatorial algorithms for shortest-path weights. Sections [TTVI12I show how 
to find unweighted /-factors in general graphs. Section [13] shows how to compute the weights of 
perturbations of /-factors. Finally Section [T4l concludes by posing several new open problems. 

Given the length of this submission, we remark that various portions can be read independently. 
Section [4] gives the whole development in the bipartite case. Section [5] gives the combinatoric 
part of the 6-matching algorithm starting from first principles. Section [6] derives the shortest-path 
structure from first principles, entirely in the language of shortest paths rather than matching. The 
related material in Section [9] is itself independent, given the definition of the shortest-path structure 
(Section 16. 1|) . Alternatively Sections [SHS] can be skipped to go directly to the combinatoric part of 
the /-factor algorithm (Sections [^^8]). 

2 Problem definitions 

The symmetric difference of sets is denoted by ©, i.e., A (B B = (A — B) U (B — A). We use a 
common convention to sum function values: If / is a real- valued function on elements and S is a 
set of such elements, f{S) denotes ^{f{v) : v G S}. Similarly if z is a function on sets of elements 
then z{S -.S € S} denotes Y^{z{S) -.SeS}. 

Let G = {V,E) be an undirected graph, with vertex set V = {1,... ,n}. We sometimes write 
V{G) or E{G) to denote vertices or edges of graph G. A walk is a sequence A = vq, ei, -yi, ... , ek,Vk 
for vertices Vi and edges ej = Vi-iVi. The notation v^v^-walk provides the two endpoints. The 
length of A is /c, and the parity of k makes A even or odd. A trail is an edge-simple walk. A circuit 
is a trail that starts and ends at the same vertex. The vertex-simple analogs are path and cycle. 

In an undirected multigraph G = (V^E) each edge e £ E has a positive multiplicity fJ.{e). Each 
copy of a fixed e £ E is a distinct edge, e.g., a trail may up to //(e) distinct copies of e. For a 
multiset S, 2S denotes S with every multiplicity doubled. Similarly 2G denotes the multigraph 
{V, 2E). If every multiplicity of S is even then S/2 denotes S with every multiplicity halved. 

For a set of vertices S C 1/ and a subgraph H of G, 5{S, H) (7(6", H)) denotes the set of edges 
with exactly one (respectively two) endpoints in S (loops are in 7 but not 5). d{v,H) denotes the 



degree of vertex v in H. When referring to the given graph G we often omit the last argument and 
write, e.g., 6{S). 

An edge weight function w assigns a numeric weight to each edge. For complexity bounds we 
assume the range of w is [— VF..Vl^], i.e., the set of integers of magnitude < W. The weight of edge 
set F (^ E is w{F) = X^eeF^(^)- ^ ^® conservative if there are no negative weight cycles. For 
multigraphs we sometimes write w{e, k) to denote the weight of the kth copy of edge e. 

Let G = {V, E) be a multigraph. For a function f : V ^ Z_|_, an f -factor is a subset of edges 
F '^ E such that d(v, F) = f(v) for every v £ V. Let G = {V, E) be a graph, where E may contain 
loops vv but no parallel edges. For a function b : V ^ Z^ , a (perfect) b-matching is a function 
X : E ^ Z_|_ such that X^^.jj^g^; x{vw) = b{v). A maximum f -factor is an /-factor F with maximum 
weight w{F). Similarly a maximum b-matching is a perfect 6-matching of maximum weight. 

To simplify the time bounds we assume matrix multiplication time is J7(n^logn). This allows 
us to include terms like 0{m) and 0(m log n) within our overall bound 0{(f)'^). Observe that 
m = 0{n'^) = 0{(jP') if there are no parallel edges, or if all copies of an edge have the same weight. 
In the most general case - arbitrary parallel edges - we can assume m = 0{n(l)) = 0{(j)'^) after 
linear-time preprocessing. In proof, for each edge uv the preprocessing discards all but the f{u) 
largest copies. This leaves < '^{d{u)f{u) : u G V} < ^{nf{u) : u G V} < n(f) edges in G. 

3 Algebraic preliminaries 

One of the fundamental ideas of our algorithms is to encode graph problems in matrices, in such 
a way that determinant of a matrix is (symbolically) non-zero if and only if the problem has a 
solution. The Schwartz-Zippel Lemma [351 [28] provides us with an efficient non-zero test for such 
symbolic determinants. For our purposes the following simplified version of it suffices. 

Corollary 1 (Schwartz-Zippel). For any prime p, if a (non-zero) multivariate polynomial of degree 
d over Zp is evaluated at a random point, the probability of false zero is < d/p. 

In order to use this lemma, we will choose primes p of size ©(n'^) for some constant c. We note 
that in a RAM machine with word size G(logn), arithmetic modulo p can be realized in constant 
time. 

However, finding the right encoding is just the first and the easiest step, whereas the more 
complicated part is to extract the actual solution from this encoding. In order to do it we use the 
following algebraic tools. The first one allows us to update the inverse of the matrix after we have 
changed the matrix itself. 

Lemma 2. (Sherman-Morrison- WOODBURY Formula). Let A be n x n non-singular matrix, 
and U,V be n x k matrices, then 

• A -\- UV"^ is non-singular if and only if the k x k matrix Ij^ -\- V'^ A~^U is non-singular and 
dei{A + UV"^) = det(A) det(4 + V'^A-^U), 

• if A + UV'^ is non-singular then {A + UV'^)~'^ = A^^ - A^^U{Ik + V^A^^U)-^V^A~K 

When k = 1 the matrices U and V become length n vectors. Such an update is called rank-one 
update. In this special case the above lemma is called Sherman-Morrison formula. Observe that for 
A; = 1 we can compute {A -|- UV'^)~^ from A~^ in 0{n^) arithmetic operations. 

In our algorithms we use the above formula but restricted to submatrices. Let R (C) denote 
set of rows (respectively columns) of matrix A. We denote the submatrix of A restricted to rows 

7 



in R and columns in C by j4[ii, C]. Harvey \\.7\ Corollary 2.1] has observed that ^"^[S, S] can be 
computed in 0(|S'|'^) time after updates have been made to submatrix A[S,S]. 

The final tool is rather recent result in symbolic computation by Storjohann |31| . He has shown 
how to compute a determinant of a polynomial matrix, as well as, how to solve a rational system 
for polynomial matrix. 

Theorem 3 (Storjohann '03). Let K he an arbitrary field, A € K[yY^^^ a polynomial matrix of 
degree W , and b € i^[y]"'^^ a polynomial vector of the same degree. Then 

• rational system solution A~^b (Algorithm 5 fSl^). 

• determinant det{A) (Algorithm 12 \31^). 

can be computed in 0{Wn^) operations in K , w.h.p. 

4 Outline of the paper, and the bipartite case 

This section has two purposes. First it presents our algorithms for bipartite graphs, a simplification 
of the general approach. (But even in the bipartite case our techniques were not previously known.) 
Second, it gives a guide to the entire paper: Each time we introduce a construction we discuss 
how it can be extended to the general case. The section ends by giving our algorithms for vertex- 
capacitated flow. 

4.1 Determinant formulations 

Consider a simple bipartite graph G, with both vertex sets Vo,Vi numbered from 1 to n. Let 
(j) = '^i fii)/'^- Define a (/) x (/) matrix B{G), the symbolic adjacency matrix ofG,f, as follows. 
A vertex i £ Vq is associated with f{i) rows, which are indexed by a pair i,r, for < r < f{i). 
Similarly j E Vi is associated with f{j) columns indexed by j,c, for < c < f{j). B{G) uses 
indeterminates x/' and yl^ and is defined by 

fi(r\ _i^ryc^ ijeE, 

B{G)i^r,j,c - Sf. ., ■ (1) 

I U otherwise. 

Observe that each edge in the graph is represented by a rank-one submatrix given by the product of 
two vectors x'^^y'^^) . Before we prove the main theorems we make the following observation that 
will show each edge can be used only once. 

Lemma 4. Let A he a symbolic n x n matrix, let R be the set of m rows, and G be the set of m 
columns of A. If A[R,G] has rank bounded hy r then each term in the expansion ofdet{A) contains 
at most r elements from A[R, G]. 

Proof. Using Laplace expansion we expand det(^) into m x m minors that contain all m rows of 
R, i.e., 

det{A) = Yl sgn(M) det{A[R, M]) det{A[Vo -R,Vi- M]), 

MCVi,\M\=m 

where sgn(M) = HceA/l"-'-)'^- Consider now each element of the above sum separately. If M 
contains > r columns of G then det(^[i?, M]) = 0, so the elements contributing to det(^) have < 
columns of G. Moreover, A[Vo - R,Vi - M] has no rows of A[R,G], so det{A[R, M]) det{A[Vo - 
R, Vi - M]) has < r entries from A[R, G]. ■ 



The determinant of B{G) is the sum of many different terms each containing exactly (j) occur- 
rences of variable pairs Xryi- Each pair Xryi corresponds to an edge ij S E. For a term a let ¥„ 
denote the multiset of edges that correspond to the variable pairs in a. Define J- to be the function 
that maps each term a to Ff,. 

Theorem 5. Let G he a simple bipartite graph. The function T from terms in det(i?(G)) is a 
surjection onto the f -factors of G. Consequently, G has an f -factor if and only ifdet{B{G)) 7^ 0|j 

Proof. First we show that the image of J-" contains all /-factors of G. Suppose F is an /-factor in 
G. Order the edges of -F that are incident to each vertex arbitrarily. If ij £ F is the r + 1st edge 
at i and the c + 1st edge at j then it corresponds to entry B{G)i^r,j,c- Thus F corresponds to a 
nonzero term a in the expansion of det{B{G)). Observe that entries B{G)i^r,j,c define a in a unique 
way, so no other term has exactly the same indeterminates. Of course there can be many terms 
representing F. 

Can det(i?(G)) contain terms that do not correspond to /-factors? We show the answer is no. 
Suppose det{B{G)) 7^ and take any term a in the expansion of det(-B(G)). a corresponds to an 
/-factor unless more than one entry corresponds to the same edge of G. This is impossible because 
edges are represented by rank-one submatrices and Lemma H] shows elements of such submatrices 
appear at most once. ■ 

Now let G be a bipartite multigraph. Let fi{e) denote the multiplicity of any edge e G E. 
Redefine the corresponding entry in B{G) by 

M(«i) 

k=l 

In other words, now edge ij of multiplicity fi{ij) is represented by a submatrix of rank fJ.{ij). Hence 
Lemma [4] shows edge ij can appear in a term of det(-B(G)) at most fi{ij) times. This leads to the 
following generalization of Theorem [5j 

Corollary 6. Let G be a bipartite multigraph. The function T from terms in det(i?(G)) is a 
surjection onto the f -factors in G. Consequently, G has an f -factor iff det{B{G)) ^ 0. 

The final point of this section concerns the complexity of using the B{G) matrix. As in most 
algebraic algorithms we evaluate B{G) using a random value for each indeterminate to get a matrix 
B. If G is a simple graph this is easily done in time 0{(j)^). But the situation is less clear for 
multigraphs. The most extreme case is exemplified by 6-matching. Consider an arbitrary edge 
ij and let B[I^J] denote the f{i) x /(j) submatrix of B that represents it. ij has multiplicity 
/i(ij) < min{/(i), /(j)}. When f[i) and f{j) are 0((/>) the time to compute the expression of ([2]) 
is O(^). So computing the QiyCJp') entries of i?[/, J] uses time B((/)^). But this can be avoided and 
we can construct B in time 0(0^), as follows. 

As observed above B[L, J] is the product of an f{i) x /u(ij) matrix X of indeterminates x^r' and 
an /i(ij) X /(j) matrix Y of indeterminates y^c' ■ Wlog assume f{i) < f{j). We have the trivial 
bound fj,{ij) < f{i). Break up the product XY into products of ^{ij) x ^{ij) matrices, i.e., break 



^The second part of the theorem suffices for undirected /-factors. But we will need the stronger claim of surjectivity 
for the weighted case. 



X into f{i)/fj.{ij) matrices of size fi{ij) x fi{ij), and similarly Y, to get XY as a sum of ^77^(7^ 
products of /i(ij) x //(ij) matrices. Using fast matix multiplication on these products, the total 
time to compute B is bounded by a constant times 

f{i)<fij) /(i)</(i) ^ /(i)</(i) » 

Generalizations The generalization of these ideas to non-bipartite graphs is given in Section [TT] 
We combine the above idea, submatrices with bounded rank, with the idea of Tutte to construct 
a skew-symmetric matrix. Additionally we need to take care of self-loops in multigraphs. In the 
Tutte matrix self-loops do not need appear since they cannot be used in a 1-factor. 

4.2 Finding /-factors 

This section gives our algorithm to find an /-factor of a bipartite multigraph. We follow the devel- 
opment from |21) : We start with an 0(i;^'^)-time algorithm. Then we show it can be implemented 
in 0{4>^) time using the Gaussian elimination algorithm of Bunch and Hopcroft [3|. 

An allowed edge is an edge belonging to some /-factor. For perfect matchings the notion of 
allowed edge is easily expressible using the inverse of B{G): ij is allowed if and only if B{G)^, is 
non-zero |24| . We will prove a similar statement for bipartite /-factors. (But such a statement fails 
for non-bipartite graphs - see Appendix lAl) For a given / define fij to be 

\f{v) - 1 ■ifv = iorv=j, 
\f{v) otherwise. 

Lemma 7. Let G be a bipartite multigraph having an f -factor. Edge ij £ E is allowed if and only 
if G has an fij- factor. 

Proof. The "only if" direction is clear: If F is an /-factor containing ij, then F — ij is an /jj-factor. 

Conversely, suppose E does not contain the chosen edge ij. Take an /jj-factor E' that maximizes 
l-F' n -F|. E' (B E contains an alternating ij-trail P that starts and ends with edges of E. In fact P 
is a path. (Any cycle C in P has even length and so is alternating. This makes E' (BG an /jj-factor 
containing more edges of E than E' , impossible.) 

ij is not the first edge of P {ij ^ E). So ij ^ P, since P is vertex simple. Thus {E (B P) + ij is 
an /-factor containing ij. ■ 

Lemma 8. Let G be a bipartite multigraph having an f -factor, let B{G) be its symbolic adjacency 
matrix, and let i £ Vq and j G Vi. Then (i?(G)~^)j^o,i,o ¥" if and only if G has fij-factor. 

Proof. Observe that 

-l)"(j'0)+"(^'0)det(S(G)^'°'J'°) 



iBiGr%^ 



0,i,0 



det{B{G)) 



where B{Gy'^'^'^ is the matrix B{G) with i,0'th row and j, O'th column removed, and n{i,k) is the 
actual index of the row or column given by the pair i,k. We have det{B{G)) 7^ since G has an 
/-factor. Hence (5(G)-i)j,o,i,o / if and only if det(5(G)*'°'J'°) / 0. Furthermore B{Gy'°'^'° is 
the symbolic adjacency matrix obtained from G for /j j-factors. ■ 
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Observe that by the symmetry of the matrix B(G), when (i?(G)~^)j^o,i,o 7^ then as well 
(-B(G)~^)j.K,i,(, 7^ for all < t < /(i) and < k < f{j)- Combining the above two lemmas with 
this observation we obtain the following. 

Corollary 9. Let G be a bipartite multigraph having an f -factor, and let B{G) be its symbolic 
adjacency matrix. The edge ij ^ E is allowed if and only if {B{G)~^)j^K,i,L 7^ for any < l < f{i) 
and < K < f{j)- 

Being equipped with a tool for finding allowed edges we can now use the Gaussian elimination 
framework from |21) . The following observation is useful. 

Lemma 10 (|21|). Let A be a non-singular (j)X (/) matrix and let 1 < i,j < n be such that {A~^)j^i ^ 
0. Let A' be the matrix obtained from A~^ by eliminating row j and column i using Gaussian 
elimination. Then A' = {A^'^)~^ (i.e., A' is the Schur complement of {A'^y'"^). 

The above lemma can be used to obtain the following algorithm that finds an /-factor. 

Algorithm 1 An 0((/>'^) time algorithm to find an /-factor in a bipartite multigraph G. 

1: 
2: 
3: 
4: 
5: 
6: 
7: 



9 
10 
11 
12 
13 
14 
15 
16 



Let B(G) be the (p x (p matrix representing G, f 

Replace the variables in B{G) by random values from Zp for prime p = 0((/)^) to obtain B 

If B is singular return "no /-factor" 

(with probability >1—-t matrix B is non-singular when B{G) is non-singular) t> by Cor. [T] 

F:=0 

Compute N := B-^ 

(each column i, l of B^l has an allowed edge, since BB^^ = I gives j, k with Bi^i^j^^Bj^ i t 7^ 0) 

for i = [l..n] do 

for i = [0..f{i) - 1] do 

Find j, K such that ij € E — F and Nj^f^^i^^ 7^ l> by Corollary [9] edge ij is allowed 

Eliminate the j, K'th row and the i, i'th column from N \> using Gaussian elimination 
(Lemma 1101 shows N = (^B^'^'^''^)~^, i.e., N encodes /jj-factors, but see below) 
Set F := F-^ ij 
end for 

end for 

Return F 



The comment of line 12 is adequate for l = 0. However l > requires an additional observation. 
To see this first recall the logic of each iteration: Let /' be the residual degree requirement function, 
i.e., the current F, enlarged with an /'-factor of the current graph, gives an /-factor of G. In line 
10, the /'-factor F' that contains ij is a subgraph of the graph corresponding to (the current) N 
and its corresponding matrix B. Now suppose the iteration for t = adds edge ip to F. When 
the row and column for ip are deleted from B, the remaining rows for vertex i still contain entries 
corresponding to edge ip (recall the definition of B{G)). So when the iteration for t = 1 chooses its 
edge ij, the corresponding /'-factor F' may contain edge ip. But F cannot be enlarged with F' , 
since that introduces two copies of ip. The same restriction applies to iterations for l > 1, but now 
it concerns all previously chosen edges ip. 

Actually there is no problem because we can guarantee an /'-factor avoiding all the previous 
ip^s exists. The guarantee is given by the following corollary to Lemma [7l (Note when r = the 
corollary is a special case of the lemma. Also, the converse of the corollary holds trivially.) 
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Corollary 11. Consider a set of edges P = {ipi, ... ,ipr}, t > 0. Suppose G — P has an f -factor. 
If G has an fij-f actor for some edge ij ^ P then G — P has an f -factor containing ij. 

Proof. The proof of Lemma [7] applies, assuming we start by taking F to be the assumed /-factor. 



Finahy observe that Algorithm 1 is implementing Gaussian elimination on B^^, the only dif- 
ference being that pivot elements are chosen to correspond to edges of the graph. If there exists 
an /-factor, there is an allowed edge incident to each vertex. Hence, even with this additional 
requirement Gaussian elimination is able to find a non-zero element in each row of B~^. 

Bunch and Hopcroft [^ show how to speed up the running time of Gaussian elimination from 
0{(j)^) to 0{(j)^), by using lazy updates to the matrix. Let us divide the columns of the matrix into 
two almost equal parts. Let L denote the first \(j)/2\ columns that are to be eliminated, whereas 
let R denote the remaining [(p/2\ columns. Bunch and Hopcroft observed that columns in R are 
not used until we eliminate all columns from LO Hence all updates to columns in R resulting from 
elimination of columns in L can be done once using fast matrix multiplication in Odi?!"^) time. By 
applying this scheme recursively one obtains an 0{<j)^) time algorithm. 

Generalizations Section [12] gives 0(</>'^) time algorithms for finding /-factors in non-bipartite 
multigraphs. There are several things that need to be done differently. As discussed above the 
criteria for finding allowed edges - Corollary [9]- does not work any more. We need to work with the 
weaker notion of removable edges, i.e., the edges that can be removed from the graph so that it still 
contains an /-factor. This forces us to use a different approach, based on Harvey |17| . which works 
with removable edges. This poses a new challenge, to handle multiple copies of edges in multi- 
graphs, as removing such edges one by one could require 0{4>^) time. To overcome this problem we 
use binary search with the Sherman-Morrison- Woodbury formula to remove multiple copies of an 
edge in one shot. 

4.3 Weighted /-factors 

In this section we discuss how to find a maximum /-factor in a weighted bipartite graph. For the sake 
of simplicity we assume in this section that the weight function is non- negative ,i.e., w : E ^ [0..W]. 
If this is not the case we can redefine w{ij) := w{ij)-\-W , what changes the weight of each /-factor by 
exactly Wf{V)/2 = W(j). Let us start by recalling the dual problem for maximum /-factors. In this 
problem each vertex v is assigned a real- valued weight y{v). We say that the dual y dominates the 
edge uv £ E when y{u) -\-y(v) > w(e), or it underrates the edge uv £ E when y{u) +y{v) < w{uv). 
The objective that we need to minimize in dual problem is 

y{V,E) = J2f iv)y (.'")+ ^ w{uv)-y{u)-y{v). 

v(^V uvdE is underrated 

The dual y minimizes yiV, E), when there exists an /-factor F such that F contains only underrated 
edges, whereas its complement contains only dominated edges. Observe that when we are given the 
minimum dual y, then the above /-factor F is a maximum weight /-factor. On the other hand, in 



®In their paper the elimination proceeds row by row, whereas it is nowadays more usual to present Gaussian 
elimination on columns. 
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order to construct such maximum /-factor we need to take into it every strictly underrated edge 
and arbitrary tight edges, i.e., edges uv £ E for which y{u) +y{v) = w{uw). Hence, we can observe 
the following. 

Lemma 12. Given an optimal dual function y, a maximum f -factor of a bipartite multigraph can 
be constructed in 0{(f>^) time. 

Proof. Let U be equal to the set of underrated edges with respect to y. Set f'{v) = f{v) — d{v, U). 
Using Algorithm [1] find an /'-factor T over the set of tight edges with respect to y. The maximum 
/-factor is equal to multiset sum of U and T, i.e., to t/ tt) T. ■ 

This lemma shows that given an algorithm for finding unweighted /-factors all we need to know 
is an optimal dual. Such an optimal dual can be obtained from the combinatorial interpretation 
as given in |14) . Let us define G^ to be G with additional vertex s G Vi and new weight edges 
sn, for all u G Vq. In G"*" we set /(s) = 1. Let /^ be the degree constraint function defined to be 
identical to / except for fyiv) = f{v) + (—1)*, where t; G V^. Let F^ be a maximum /^,-factor in 
G"*". To show that F^ always exists take F to be any /-factor in G. Now, when -y G Vq then F -\- sv 
is an /^-factor, whereas when v ^Vi — s then for any uv (z F, F — uv + su is an /^-factor 

Theorem 13 ([H]). For a bipartite multigraph with an f -factor, optimal duals are given by y{v) = 
{-iyvu{F,)forveVi. 

Hence, in order to construct optimal dual we need to know weights w{Fi,), for all t; G Vq U Vi. 
At first sight it might seem that we did not gain anything, as instead of finding one F factor now 
we need to find all factors F^,. However, we do not need to find these factors. We only need to know 
their weights, which is much easier. And the following lemma shows that we just need to know 
w{Fi,) for one side of the bipartite graph. 

Lemma 14. For a bipartite multigraph with an f -factor, let y{v) be an optimal dual for each 
V €z Vi. An optimal dual y{u) for u ^ Vq is equal to the largest value y^ that makes at least f{u) 
edges incident to u underrated, i.e., \{uv^^ : yu < w{uv) — y{v)^\ > f{u). 

Proof. Observe that there are at least f{u) underrated edges incident to u with respect to optimal 
dual y, as each maximum /-factor needs to contain only underrated edges. On the other hand, the 
complement contains at least d{u) — f{u) dominated edges. This fixes the largest possible value for 
y{u) as the value given in the lemma. ■ 

Now consider a simple bipartite graph G and similarly to ([T|) define B{G) as 

^[^)i,r,j,c — N 

I otherwise, 

where 2 is a new indeterminate. Theorem [5] shows that there is a mapping F from terms of det{B{G)) 
onto /-factors in G. Consider a term a in det{B{G)). Observe that its degree in z is equal to the 
weight of F{a) because the powers of z get added in the multiplication. For a polynomial p, denote 
the degree of p in z by deg^(p). We obtain the following observation. 

Corollary 15. For a simple bipartite graph G, deg^(det(-B(G))) equals the weight of a maximum 
f -factor in G. 
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To compute /^-factors in G^ we use the following auxiliary graph. Let G* be G"*" with an 
additional vertex t G Vq that is joined to every vertex oi Vi — s by an edge of weight zero. Set 
fit) = 1. 

Lemma 16. Every v ^Vi satisfies deg^{a,di]{B{G^))i,fl^tfl) = w{Fy). 

Proof. Observe that 

adj(i?(G,)),,o,i,o = (-l)"(*'0)+"(-'0) det(i?(G,)*'°''''°), 

where i?(G,,)*''^'''''' is the matrix B{G<t:) with row t,0 and column v,0 removed and n{i,r) gives the 
order of rows and columns indexed by pairs i, r. By Theorem[5]we know that det(i?(G*)) consists of 
terms corresponding to /-factors in G*. Hence the above equality shows terms of adj(B(G^,))t,^o,t,o 
correspond to /-factors that use edge tv, but with this edge removed. These are exactly the /„- 
factors in G+, because forcing the /-factor to use edge tv effectively decreases f{v) by 1. As we 
observed in Corollary [15] the degree of z equals the total weight of corresponding /-factor. The 
lemma follows. ■ 

Recall that the adjoint of a nonsingular matrix A is dei{A)A~^ . The lemma shows we are inter- 
ested in column t, of the adjoint. So let ej^o be a unit vector, with 1 in row t, and zeroes elsewhere. 
Then the desired weights are found in the vector adj(-B(G^,))ef^o = det(i?(G*))i?(G*)~^et^o- This 
leads to the following algorithm to find optimal duals for weighted bipartite /-factors. 

Algorithm 2 An 0{W(j)''^) time algorithm to find optimal duals y in a bipartite graph G, /. 
1: Let B[G^) be x (p matrix representing G* 

2: Replace x and y variables in B{G^:) by random values from Zp for prime p = 0(0'^) to obtain B 
3: Compute vector a := adj(i?)et^o = det(i?)i?~^et^o- ■> requires OiWcj)'^) time using Theorem [3] 
4: for u € Fi do 

5: (deg^(a„) = deg^(adj(S(G*))^,o,t,o) with probability > 1 - ^) > by Cor. [U 

6: Set w{Fy) := deg^(a^) > equality holds by Lemma [TBI 

7: Set y{v) := w{Fy) > y{v) is optimal by Theorem [T3] 

8: end for > by union bound all y{v) are correct with probability > 1 — ^ 

9: for u E Vb do 

10: Set y{u) := maxjy^ : \{uv^'^ : yu < w{uv) - y{v)}\ > f{u)} 

11: > y{u) is optimal by Lemma [Ml 

12: end for 

Combining the above algorithm with Lemma [12] we obtain an 0{Wn'^) time algorithm for 
maximum /-factors in weighted bipartite graphs. It can be observed that this development works 
for bipartite multi-graphs as well, when one changes the definition of B[G) to 

M(y) 

k=\ 

As in Section [31 some problem instances require a more careful construction of i?. But the algorithm 
of Section [4] is easily extended. First note that 6-matching problems can be handled using the 
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algorithm of Section U] unchanged, since every copy of a fixed edge ij has the same weight w{ij, k) = 
w{ij). The most general case (/-factor problems with high multiplicities and parallel edges of 
different weights) is easily handled as follows. Decompose G into multigraphs G^ containing the 
edges of weight w, for w = 0, ... , W. So B{G) = Ylw=o z'"B{G^), and analogously, B = ^^^^ z'^B^ 
(where all -B^'s use the same random values for the x and y variables). Compute each B^ using the 
algorithm of Section [4] and combine. Each Byj is found in time 0{(f)^) so the total time is 0{W(j)'^). 

Generalizations The relation between the primal and dual problems in the bipartite case is 
considerably simpler than the general case. The latter not only contains dual variables for vertices, 
but also for subsets of vertices. Such subsets with non-zero dual value are called blossoms. We 
can prove these blossoms are nested and so form a blossom tree. Moreover, for each blossom we 
need to find a set of spanned edges and incident edges that are all underrated. This cannot be 
done in such a simple way as Lemma [TH But again knowing the weights of maximum /^-factors 
turns out to be enough. The procedure that deduces all this information is highly nontrivial and 
is described in Section [HI This section is preceded with two special, simpler cases. First Section [5] 
considers maximum 6-matchings, where there are no underrated incident edges. Second Section [6] 
considers shortest paths in undirected graphs with negative weights, where each blossom has exactly 
one underrated incident edge. Finally Section [T3] generalizes Algorithm [2] to compute w{Fi,) for 
non-bipartite graphs. Alternatively for the special case of shortest paths w{F^) can be computed 
combinatorially as in Section [lOl 

4.4 Min-cost max-flow 

This section presents network flow algorithms. We discuss maximum flows and minimum cost 
maximum flows, both in vertex-capacitated networks. 

We are given a directed network N = {V, E), with source s and sink t, s,t G V. For convenience 
let V~ denote the set of nonterminals, V — {s,t}. The edges and nonterminal vertices have integral 
capacities given by c : V^ UE — > [1..D]. Let g : V xV ^ Tjhe a flow function. Besides the standard 
edge capacity and flow conservation constraints we have vertex capacity constraints, i.e., for each 
vertex v ^ s,t we require 

J29{u,v) < c{v). 
uev 

We begin by constructing a bipartite graph Gn whose maximum /-factor has weight equal to 
the value of a maximum flow in N. Wlog assume that no edge enters s or leaves t. The construction 
proceeds as follows: 

• for each v G V~ place vertices Vin, Vout in Gn] 

• also place vertices Sout, tin in Gn; 

• for each v S V~ add c{v) copies of edge VinVout to Gn; 

• for each {u,v) G E add c{u,v) copies of edge UoutVin to Gn', 

• ioxv £V~, set f{vin) = f{vout) = c{v); 

• set f{Sout) = fiUn) = c(y~)] 

• set w{e) = 1 for each edge e leaving Sout and w{e) = for every other edge of Gn; 

• add c{V~) copies of edge SoutUn to Gjy, all of weight 0. 

Note that in addition to the weight edges Soutiin, Gn may contain edges Souttin of weight 1 
corresponding to an edge st G E. 
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Figure 1: A vertex v of capacity c{v) = 3 in A^ is represented in G^ by two vertices Vin and Vout 
connected by 3 edges. The /-factor in (b) is marked with soUd edges. Observe that the /-factor 
must choose the same number of edges going in and out of v. 
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Corollary 17 (|15|). Let N be the flow network. The weight of the maximum f -factor in Gn is 
equal to the maximum flow value in N. 

Proof. The main idea of the reduction is shown in Figure [1] Observe that /-factors in Gjy correspond 
to integral flows in N that fulfill the flow conservation constraints. Moreover the edge capacities are 
not exceeded since an edge cannot be used by an /-factor more times than its capacity. Similarly 
a vertex cannot be used more times than its capacity. The only edges with non-zero weights are 
edges incident to Sout, so the maximum /-factor maximizes the amount of flow leaving s. Finally 
all this flow must wind up at t, since any v ^ s,t has an equal number of /-edges incident to Vin 
and Vout- ■ 

Observe that Gn has f(y) < 4:c{V~). So the algorithm of Section |13] uses d{c(y~)'^) = 
0{{Dn)^) time to find a maximum flow in a network. 

Now assume that an edge (u, v) of N has a cost a^^^, G [— M^..VF], i.e., the cost of sending t units 
of flow on edge (n, v) is linear and equals tOu^v First flnd the maximum flow value fmax in A^. Then 
modify the construction of Gjy to Gjy^a in the following way: 

• add vertices Sj„ and tout] 

• add c{V~) copies of edges SinSout and tintout; 

• set f{Sin) = f{tout) = C{V~) and f{Sout) = fiUn) = c(y~) + fmax] 

• for each copy of the edge UoutVin set w{uoutVin) = -au,w 

In Corollary [T7] we observed that /-factors in Gn correspond to feasible flows in A^. The new 
vertices Sj„ and tout allow an /-factor to model flow along cycles containing s and t. The new values 
of / at the terminals make the corresponding flow from s to i have value fmax- Thus /-factors 
in Giv.a correspond to maximum flows in A^, and since costs of edges are negated between both 
networks, we have the following. 

Corollary 18 (|15|). Let N he a flow network with linear edge costs. A maximum f -factor in GN,a 
has weight equal to the minimum cost of a maximum flow in N. 

Now the algorithm of Section [4.31 gives an 0{Wc{V~)^) = 0{W{Dn)''^) time algorithm to find 
a min-cost max-flow. 

Let us extend this reduction to convex cost functions. Assume the cost of sending t units 
of flow on an edge {u,v) is given by a convex function au^v{t) such that marginal costs satisfy 
Tnu,v{t) = 0'u,v{t) — OLu,v{t — 1) G [— M^..VF]. As usual for this scenario [T] we assume au^v is linear 
between successive integers. This ensures that there exists an integral optimal solution. We encode 
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such cost functions in the graph by assigning different costs to each copy of an edge, i.e., the kth. 
copy of edge UoutVin has cost w{uoutVin,k) = —mu^v{k) for 'i<k< c{u,v). 

Lemma 19. Let N be a flow network with convex edge costs. A maximum f -factor in G^.a has 
weight equal to the minimum cost of a maximum flow in N . 

Proof. The maximum /-factor in G^^a, when using t copies of edges between Uout and Vin, will 
use the most expensive copies. The convexity of a^^y implies mu,v{t) is a non-decreasing function. 
Hence we can assume that edges with costs — m,„^t,(l), ... , — mu,j;(i) are used. These costs sum to 
-au,v{t)- ■ 

Clearly this extension does not change the running time of our algorithm. So we find a min-cost 
max-flow for convex edge costs in time 0{W{Dn)'^). 

It is also easy to incorporate lower bounds on flow into the reduction. Suppose that in addition 
to the upper bound function c we have a lower bound function I : V~ U E ^ [0..D]; (. restricts the 
flow on edges and flow through nonterminal vertices in the obvious way. To model £, change the 
multplicity of VinVout to c{v) — i{v) and the multiplicity of UoutVin to c{u, v) — i{u, v); in addition for 
V £ V~, f{vin) becomes c{v) — X^„^('W, ?;) and f{vout) becomes c{v) — X]tj^(f,ii). It is easy to see 
that the correspondence between /-factors and feasible flows is maintained, as is the correspondence 
between weights and costs. (In particular starting from an /-factor, adding i{u, v) copies of each 
edge UoutVin gives a subgraph that obviously corresponds to a flow satisfying all lower bounds on 
edges. It is easy to check this flow also satisfies all lower bounds on vertices.) 

The construction can also be generalized to bidirected flows. In directed graphs an undirected 
edge can have two orientations, in-out and out-in, whereas a bidirected graph allows four possible 
orientations in-in, out-in, in-out and out-out. (We also allow loops, especially of type in-in and 
out-out.) A bidirected network G gives a non-bipartite graph Gi\f. However as we show in the 
remainder of this paper, the /-factor problem for non-bipartite graphs can be solved in the same 
time bounds as the bipartite case. 

5 Weighted 6-niatching 

This section gives our algorithm to find a maximum 6-matching. b : V ^ Z-^ can be an arbitrary 
function. But we remark that our algorithm is of interest even for the case 6=1 (ordinary matching): 
It achieves the same time bound as |7], and is arguably simpler in both derivation and algorithmic 
details jj 

We start with two remarks that modify the definition of the problem. First it is most often 
convenient to use the language of multigraphs. In this view we think of G as a multigraph with an 
unlimited number of copies of each edge. Specifically each e £ E has 1 + max^ b{v) copies, each 
with the same weight w{e). (Also recall from Section [2] that G may have loops.) A 6-matching is a 
subgraph of this multigraph. 

Second, our algorithm actually works on critical graphs, defined as follows. Given a function 

6 : y — )• 1i~^, for each v £ V define b^ : V ^ Z_|_ by decreasing b{v) by 1 (keep all other values 
unchanged). A graph is b-critical if it has a perfect 6j,-matching for every v £ V. Given a 6-critical 



^The reader may enjoy working through the derivation for the matching case b = 1. There will be no need to 
work with a multigraph. This endeavor will show how the two problems differ, especially in some details introduced 
by multigraphs. 
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graph, our algorithm produces a "blossom tree" from which, for any v £V, a maximum 6^-matching 
can be easily extracted. 

Such an algorithm can find a maximum 6-matching as follows. Suppose we seek a maximum 
6-matching on G. Assume every b{v) is positive (discard 0- valued vertices). Form G' by adding a 
vertex s with b{s) = 1, plus edges sv, v G V{G) and loop vv, all of weight 0. Any v £ V + s has 
a 5i,-matching. (For v = s take F, an arbitrary 5-matching on G. For v ^ s, for any edge vw G F 
take F — vw + ws. This is a 6„-matching even if vw is a loop.) So G' is 6-critical and a maximum 
6s-matching on G' is the desired 6-matching. 

5.1 The heaviest blossom 

This section defines the ("-value of an edge and shows how it reveals the heaviest blossom (the set 
B defined after Lemma [22]) . Consider a weighted 6-critical graph G. Recall that G may contain 
self-loops vv. All of what follows is valid even when there are such loops. Assume that each multiset 
of < ^yb{v) edges has a distinct weight, if we do not distinguish between parallel copies of the 
same edge. We can enforce this assumption by taking a very small e > and adding e* to the weight 
of every copy of the ith. edge of E. (Section 15.31 returns to the original unperturbed weights.) 

Any vertex v has a maximum 6„-matching F^. Fy is unique up to parallel copies of the same 
edge, by the perturbed weight function. Wlog assume further that any two matchings Fu and F^ 
have as many common edges as possible (i.e., for any xy £ E they use as many of the same copies 
of xy as possible). 

We start with a well-known principle. Take two vertices u, v. Call a trail alternating (for u and 
v) if as we traverse it the edges alternate between F^ — F^ and F^ — F^. 

Lemma 20. For any two vertices u,v, F^® F^ is an alternating uv -trail that starts with an edge 
of Fy — Fu and ends with an edge of Fu — Fy. 

Proof. First we show that Fu ® Fy contains an alternating trail as described in the lemma. (Then 
we show that trail constitutes all of Fu © Fy.) Let T be a maximal length alternating trail that 
starts at u with an edge of Fy — Fu- Such an edge exists since u has greater degree in Fy than Fu. 
Let T end at vertex x. 

If X ^ u,v then, since d{x,Fu) = d{x,Fy), we can extend T with an unused alternating edge. 
(This is true even if T has a previous occurrence of x.) li x = u then we can extend T with an 
unused alternating edge - the argument is the same using the facts that d(u, Fy) = d{u, Fu) + l and 
the first edge of T is in Fy — Fu- Suppose x = v. If the last edge of T is in Fy — Fu we can extend 
T with an unused edge of Fu — Fy since d{v,Fu) = d{v,Fy) + 1. The remaining possibility is that 
T ends at v with an edge of Fu — Fy, thus giving the desired trail. 

Now we show there are no other edges. Fy is the disjoint union of its edges in T and a multiset 
of edges R. (R may contain copies of edges in T.) We claim any vertex x satisfies 

d{x,R) = d{x,Fy)-d{x,TnFy) (3) 

= b{x)-\d{x,T)/2]. (4) 

Observe that this claim completes the proof of the lemma: Define the function 5' by setting b'{x) 
to the quantity of ^. Thus i? is a ^'-matching, in fact a maximum ^'-matching. By symmetry Fu 
is the disjoint union its edges in T and a maximum ^'-matching. The perturbed weight function 
implies that ^'-matching is also R. Thus Fy and Fy agree outside of T, so Fy® Fy = T. 

18 



To prove the claim, ([3]) holds by definition. For dH first observe that x ^ v implies 

d{x, F^) = h{x) and d{x, T n F^) = \d{x, T)/2] . 

In proof the first relation follows from definition of Fy. For x ^ u,v the second relation holds 
because the edges through x alternate. For x = u the second relation holds because the first edge 
of T is in Fy, and all other pairs of edges through u alternate. Substituting the two relations into 
© gives gD. 

li X = V then d{v,Fy) = h{v) — 1 and d{v,T D Fy) = (d{v,T) — l)/2. The latter holds because 
the last edge of T is not in Fy. Furthermore all other pairs through v alternate. Again substituting 
these two relations into ([3]) gives ([4]). ■ 

In an arbitrary multigraph let C be an odd circuit containing a vertex u. Choose a traversal 
of C that starts at u. Define Cu to consist of alternate edges in this traversal, omitting the first 
edge at u as well as the last. When C contains > 2 edges at u, Cu will not be unique. But suppose 
the traversal starts with edge uv, and we define Cy by the same traversal only starting at v. Then 
Cy, Cy and uv form a partition of C (since C„ starts by containing the edges of C — C^, and this 
pattern continues until the traversal reaches v). Also given any choice of a Cx, define C^ to consist 
of alternate edges of C, beginning and ending with the edge incident to x. Clearly C^ and C^ form 
a partition of the edges of C. 

To define the central concept, for any edge uv let 

C{uv) = w{Fu) + w{Fy) + w{uv). 

We shall see that ( gives the values of the optimum blossom duals as well as the structure of the 
optimum blossoms. 

Lemma 21. Any edge e of a b-critical graph belongs to an odd circuit of edges with Q-value > C{e). 

Proof. Let e = uv. Assume u ^ v else the lemma is trivial. Furthermore assume e ^ FyL) Fy. This 
is the crucial assumption! It is justified since there are b{v) + 1 copies of e; furthermore proving the 
lemma for this copy of uv proves it for every copy. 

Let T be the trail of Lemma [20l It clearly has even length. Extend it by adding a copy of 
edge e. (T may already contain a different copy of e, one in Fu ® Fy.) We get an odd circuit C. 
Traverse C by starting with the first edge of T and ending with the edge uv € C — T. We get C 
partitioned into Cu, Cy and uv. Furthermore Cu = C D Fy, Cy = C D Fy, and a set of edges R 
satisfies R = Fu — C = Fy — C. Thus 

C{uv) = w{Fy) + w{Fy) + w{uv) = w{C) + 2w{R). 

Now take any edge rs G C. Traverse C by starting with edge rs. For t G {r, s} let Ht be the 
multiset RU Ct. It is easy to see Ht is a ftj-matching by comparing it with Fy = RU Cy. Then 

C(rs) > w{Hr) + w{Hs) + w{rs) = w{C) + 2w{R). 

The two displayed equations show C is the desired circuit. (In this argument some vertices t may 
have Ct multiply defined. That's OK.) ■ 
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Let C* be the maximum value of (^. Let E* be the set of edges of (^-value (*. We shah see that 
E* is essentially the heaviest blossom and ("* its dual value. 

Lemma 22. Any edge e = vw G E* belongs to an odd circuit C '^ E* . Furthermore (i) F^, — C = 
F^-C and in) F„ D {j{v) D 6{v)) C E* . 

Proof. Lemma [21] shows the odd circuit C exists. To prove (i) recall from the proof that 

w{C) + 2w{R) = C 

and F^ = RL) Ci,, Fyj = RU Cyj . The last two equations imply (i) . 

For (ii) take any edge uv G F^. Keep C and R as already defined for vw. We can assume 
uv E R, since otherwise uv £ C^ '^ C (^ E* . Let Hu = R — uv + C" . It is easy to see H^ is a 
6„-matching by comparing it with F^ = RU C^ (note u may or may not belong to C) . Then 

C{uv) > w{Hu) + w{Fy) + w{uv) = w{C) + 2w{R) = C- 

This implies equality holds and proves (ii). ■ 

Call a connected component B of E* nontrivial if it spans at least one (^*-edge|f| {B may consist 
of a single vertex v with one or more loops vv.) The next lemma shows that B behaves like a 
blossom, i.e., it can be shrunk to a single vertex. We begin the proof with two observations. First, 
any two vertices v,w £ B have 

F„-E*{B) = F^-E*{B). (5) 

In proof, since B is connected we need only show ([5]) when vvo € E*{B). That case follows by 
applying Lemma \22[[) with C C E*{B). 
Next observe for any v £ B, 

F,n{-f{B)u6iB))CE*{B). (6) 

This follows since an edge xy in the left set but not in the right, with x £ B, belongs to Fi, — E*{B) = 
F^ - E*{B) by ©. But Lemma [22]^n) shows xy £ F^ n 5{x) C E* . This implies y G S and 
xy £ E*{B), contradiction. 

Lemma 23. For any vertex v, \Fy n 6{B)\ equals for v £ B and 1 for v ^ B. 

Proof. For v £ B this follows from ([6]). So suppose v ^ B. Choose any vertex t £ B. Lemma [20l 
shows Fy ® Ft is an alternating ut-trail. Let T be the subtrail from v to the first vertex of B, say 
vertex x, and let e be the last edge of T. e £ Fy (e ^ Ft by the lemma for t, i.e.. Ft 6{B) = 0). 
Applying ([5]) to t and x shows F^ © Fx contains T. So Lemma [20] applied to v and x shows 
Fv © Fx = T. (More precisely, the first part of the proof of Lemma 1 could have chosen the current 
T as its vx trail. The second part of the proof shows Fy © Fx has no other edges.) Thus e is the 
unique edge of Fy incident to B. ■ 



*For ordinary matching B is a cycle that comprises all of E* . 
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5.2 Iterating the construction 

The last lemma generalizes [141 Lemma 3.2] for ordinary matching. The rest of the derivation closely 
parallels ordinary matching, as described in |14) and [7]. In a nutshell, the optimum blossoms and 
their duals are found by running the shrinking procedure of |14| . The nesting of the contracted blos- 
soms gives a "blossom tree" B. Any desired maximum 6^-matching {v € V{G)) can be constructed 
by a top-down traversal of B that finds a perfect 6'-matching (for appropriate b') at each node of 
B. This last step is the biggest difference from ordinary matching: In ordinary matching the edges 
found at each node of B are simply alternate edges of a cycle. For completeness we present all of 
these remaining details, of course modified for 6-matching. 

Consider a nontrivial connected component B oi E* , as above. Contract B to a vertex B, with 
G the resulting multigraph. Assume a contraction operation can create parallel edges but not loops. 
(Thus the graph changes even when B is a. single vertex v with E* consisting of loops vv.) Extend 
the degree- const aint function 6 to G by setting b{B) = 1. Lemma [23] shows G is 6-critical. From 
now on V designates the vertex set of the original given graph. 

Next we define a weight function on G. For convenience we designate edges of G by their 
corresponding edge in G (i.e., an edge of G is written as uv where u,v £ V and possibly one of 
them belongs to B). For v G B let 

B„ = F^nj{B). 

For edge uv in G define a weight w{uv) by 

w{uv) = w{uv) u,v ^ B, 

= w{uv) + w{B^) V G B. 

This definition preserves the structure of G in the following sense. For any vertex v let v be 
its image in G. Let F-y denote the (unique) maximum 6trniatching in G. For a fixed vertex v, let 
vertex xGBhevifvGB, else the end of an edge of Fi,. x is uniquely defined by Lemma [23] 

Lemma 24. (i) F^, = F^UB^. 

(ii) Any edge of G has the same C,-value in G and G. 

Proof, (i) Let i?^ be the image of F^ in G. Lemma [23] shows Hia is a fe^j-niatching. Furthermore 
Fv = HyUBx (for V ^ B this follows from the optimality of F^). So we must show H^j = Fjj. Clearly 
it suffices to show 

w{F:u)<wiH^). (7) 

Suppose V G B. FjjU By is a 6^-matching. Since it weighs no more than F^, w{Fy) < w{Hy). 
This is equivalent to ([7]) since neither set contains an edge incident to B. 

Suppose V ^ B. The definition of w shows w{Fy) = w{Hy), as well as w{Fy) < w{Fy) (by the 
optimality of Fy). These relations combine to give ([7]). 

(a) Part (i) shows w{Fy) is W{Fy) for v ^ B and w{Fy) — w{By) for v G B. The two cases of 
(ii) (depending on whether or not the edge is in S{B)) follow. ■ 



The following shrinking procedure |14| iterates the construction of G. It also constructs a tree 
that represents the nesting of the contracted blossoms: 

Start by creating a one-node tree for each vertex of V. Then repeat the following step until 
the graph consists of one vertex: 
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Let i? be a nontrivial connected component of E* . Form a tree whose root represents B; the 
subtrees of the root are the trees that represent the vertices of B in the current graph. Then 
change the current graph by contracting B and defining b and W as described above for G. 

The construction can be iterated since as noted each G is 6-critical. Lemma [2¥ln) shows ^* 
never increases from one iteration to the next. 

CaU the final tree the blossom tree B of G, and each nonleaf a blossom of B. A child C of B 
is either a blossom,'Child or a vertex-child (i.e., a leaf of B). Note that a vertex of V can occur in 
B as a singleton blossom as well as a leaf. For any node B of B let V{B) denote the set of leaf 
descendants of B. 

The next goal is to describe how the edges of any F^, v GV, are distributed in B. We begin with 
several definitions. Consider the iteration in the shrinking procedure that ends by contracting B. 
Let Gb be the graph at the start of the iteration. Thus each child C of B corresponds to a vertex of 
B in the graph Gb- Let Ci^) be the value of (* in this iteration. Let E*{B) be the corresponding 
set of edges (i.e., the edges of (^-value C,{B) that join two vertices of B in Gb)- 

We now generalize Lemma [24lfi) to any graph Gb- Consider any vertex v € V{B). As before 
let V be the image of v in Gb (so u is a child of B). Fy denotes the unique maximum 6jr matching 
of Gb- For each child C of i? let x S V{G) be the vertex that is either v (if C = v) or is the end of 
an edge of Fy. x is uniquely defined since a blossom C has b{G) = 1, and there is no choice for x if 
C is a vertex of V. We claim 

Fyn-f{v{B)) = Fyn-f{v{B)) u U F,nj{v{C)). (8) 

C a child of B 

The import of ([8]) is that it defines the entire set F^. Specifically let V be the root node of B. 
Clearly any v £ V has F^ = F^ 7(^(V)). So the entire set Fy is defined by applying ([S]) to V and 
then recursively to the children of V. 

Also note that in ([8]) if C is a vertex-child of B then the expression Fx 'y{V{G)) is empty, by 
convention. 

To prove (J8]) first observe Fy n E{Gb) = Fy. This follows by repeated applications of Lemma 
[Mlfi). which shows the edges of Fy in G form Fy. 

The observation justifies the first term FyC] ^{V[B)) in ([8]). To complete the proof note that 
the remaining edges of Fy are contained in the various blossom-children G of B. We show these 
edges are given by the terms F^ Pi "f{V{G)) in dH). If C = U this is obvious since x = v. If C ^v 
then Fy Pi 5{G) is a unique edge incident to x. The optimality of Fy implies it agrees with Fx in 
l{V{G)). ® follows. 

dU]) shows that in the graph G^, -F^jPl ^{B) C E*{B). Applying this recursively with ([S]) shows 
Fy consists entirely of C*-edges, more precisely, in any blossom B, the edges of Fy that join children 
of B (or that are loops incident to a vertex that is a blossom-child of B) are in C{B). 

5.3 The efficient b-matching algorithm 

An efficient algorithm cannot work with the perturbed weight function. So suppose we execute the 
shrinking procedure, starting with the original unperturbed weight function w. Each iteration will 
find a connected component B that was found in the construction of i3. More precisely an iteration 
that contracts blossom G of B will be skipped iff the parent of G in B, say B, has [C(C')J = [C(^)J • 
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So the shrinking procedure wih construct a tree W that is a contraction of B (more precisely an 
edge of B from parent B to child C is contracted iff it satisfies the above relation [C(C')J = [C(-B)J)- 

For any blossom B of W, define the graph G*^ as follows. Consider the iteration of the shrinking 
procedure that contracts B. The vertices of G*^ are the vertices of the current graph that are 
contained in B. (These vertices are the children of B in W.) The edges of G^ are the edges E*{B) 
in the shrinking procedure. 

The following observation is key to our algorithm. It enables us to find the edges of any desired 
maximum 6^-matching that occur in each graph G*^. Take any vertex v S V{B). Let v be the image 
of V in G*Q. The edge set F^ n E*{B) is a fe^matching of G*g. In proof first recall that Fy consists 
entirely of C*-edges. Then apply ^ to B and to each descendant of B (in B) that gets contracted 
into B when we form W from B. 

We find a maximum 6„-matching using a recursive procedure b_match that finds the desired 
edges in each graph G*^. For a blossom B oiW and a vertex v € V{B), b_match{B,v) finds a 
6^-matching of C*-edges in the subgraph of G induced by vertices V{B). It works as follows. 

Let V be the child of B that contains v. First find a 6^r matching of G*^, say Hy. (Hjj exists 
by the above key observation.) Add Hjj to the desired set. Then complete the desired set using 
recursive calls on the children of B. Specifically for each blossom-child C oi B (in W), execute 
b_match{C, x) where x E V{C) is v (if C = v) or the end of an edge of Hy. 

As before let V be the root of W. Let Hy be the 6^-matching of the given graph G that is found 
by b _match{V , v) . In general Hy is not Fy. For instance in a recursive call b _match{G , x) , x may 
differ from the vertex x given by ([8]). But because Hy consists of ^*-edges, we can prove Hy is 
optimum using duality, as follows. 

As before for any blossom BofW let C{B) be the value C* when B is created. (Now C(-S) is an 
integer.) For any blossom i? 7^ V let p{B) be its parent in W. Let W~ be the set of all blossoms of 
W. Define functions y : V ^ Z, z : W~ — ?• Z by 

y{v) = -w{Fy) v£V, 

C(V) B = V, 



^^^^ \ C{B) - C{P{B)) B£W- -V. 

Recall our convention for summing functions (Section[2]), e.g., an edge e = uv has y{e) = y{u)-\-y{v). 
Any blossom B has Q{B) = z{A : V{B) C l/(^)}. The definition of ( shows any edge e of G has 
w{e) = y{e) + C(e)- So e G E*{B) implies 

w{e) = y{e) + C{B) = y{e) + z{A : e (^ V{A)]. (9) 

Let F be either Fy or Hy. Both sets F weigh 

Y,w{e) = Y.y{e) + z{A:e<^V{A)} 

eeF eeF 

= (6.y)(y) + j;MA):eCy(A)} 



sG-F 



biViA)) - 1 



{Ky)iV)+ Y. ''''r~' 4A). 



A€W- 



For the last line recall that for both sets F, F ^(V{A)) is a 6^-matching for some x E y{A), so 
|Fn7(y(A))| = ^ 2 • ^^ conclude w{Hy) = w{Fy), i.e., b_match constructs a maximum 
6j,-factor. 
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To show the algorithm is Las Vegas we only need to observe that the functions y, z fulfill all 
the requirements to be optimum linear programming duals. This follows from two more properties. 
First, the function z is nonnegative except perhaps on V. This follows since C* never increases. 
Second any edge e of G satisfies t(;(e) < y{e) + z{A : e C T/(A)}. To see this first note that 
C(c) < C{B) for B be the first blossom created with e C V{B). Use this relation to compute w{e) 
as in dl]). 

These two properties make y, z optimum linear programming duals. So if the properties are 
satisfied, and the algorithm finds a 6^-matching composed of edges that satisfy ([9]), that matching has 
maximum weight. (This is easily verified without appealing to linear programming, by computing 
the weight of an arbitrary 6i,-matching F similar to the computation of w{F) above.) 

To summarize the algorithm of this section works as follows. Consider a 6-critical graph G, with 
edge weights w^e). Assume we are given the weight of every maximum 6^-matching, i.e., every value 
w{F^), V £ V. Start by executing the shrinking procedure to construct the tree W. For any vertex 
f G y, to find a maximum 6„-matching call b _'match{V , v) . 

The total running time for this procedure is 0{(p'^). To prove this we will show the shrinking 
procedure uses 0{mlogn) time and h_match uses 0{(j)^) time, thus giving the desired bound. First 
note that W has < 3n nodes (n leaves, n blossoms that are singletons, and n larger blossoms). 

The shrinking procedure starts by using the given values w(Fy) to compute the (^-value of each 
edge. Then it sorts the edges on decreasing (^- value, in time 0{mlogn). 

We use a set-merging algorithm [6] to keep track of the contracted vertices. That is, for any 
V GV, find{v) gives the contracted vertex currently containing v. The operation union{A,B,C) 
merges two contracted vertices A and B into a new vertex C. 

The iteration of the shrinking procedure that creates B starts by constructing an adjacency 
structure for the C*-edges of graph Gb- (An edge uv {u,v G V) joins the vertices given by f ind(n), 
find(f).) Then it finds the nontrivial connected components. Each edge of G is in at most one 
graph Gb- So the total time is 0{m) plus the time for n unions and 0{m) finds [6]. 

b_match runs in the time to find ^jj-matchings on all the graphs G^. These graphs contain 
a total of < 3n vertices, < m edges, and degree constraints b totalling < 3(p (the vertices of V 
contribute exactly (p and each blossom contributes 1 , giving < 2n + (p < 3(j)) . So using an algorithm 
that finds a 6-factor in time 0{(p^) gives total running time 0[(f)'^). Since the algorithm is randomized 
there is a minor point of controlling the error probability when the graphs get small, this is easily 
handled. 

6 Shortest paths 

This section discusses the single-source shortest-path problem on conservative undirected graphs. 
It defines a "shortest-path tree" structure for conservative graphs and proves this structure always 
exists. The structure is a special case of the one for general /-factors. It is used as an example in 
SecEl 

6.1 The shortest-path structure 

Let (G, t, w) denote a connected undirected graph with distinguished vertex t and conservative 
edge- weight function •(/;:£'—?• M. We wish to find a shortest path from each vertex to the fixed sink 
vertex t. 
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Figure 2: Conservative undirected graph. Vertex labels are shortest-path distances; arrows show 
the first edge of shortest paths. Dashed edges are not in any shortest path. 

Bellman's inequalities needn't hold and a shortest-path tree needn't exist (e.g., the subgraph 
on {a.b.c} in Figl2]). Nonetheless the distance variables of Bellman's inequalities are optimum for a 
related set of inequalities. 

We begin by defining the analog of the shortest-path tree. When there are no negative edges 
this analog is a variant of the standard shortest-path tree (node V below). FigslJHS] illustrate the 
definition. In Figl2]the arrow from each vertex v gives the first edge in f 's shortest path. This edge 
is e(f) in the definition below. More generally this edge is e(A^) for any node A^ that it leaves, e.g., 
ce = e(c) = e({a, 6, c}). 

Definition 25. A generalized shortest-path tree (gsp-tree) T is a tree whose leaves correspond to the 
vertices ofG. For each node N ofT, V{N) denotes the set of leaf descendants of N, V{N) C V{G). 
Let V be the root of T ■ For each node N, V{N) contains a sink vertex denoted t{N); for N ^ V, 
t{N) is the end of an edge e{N) S 5{V{N)). For N = V the sink is t, and we take e(V) = 0; for 
N ^ V, t{N) and e{N) are determined by the parent of N as described below. 

Consider an interior node N ofT, with children Ni, i = 1, ... ,k, k > 2. V{Ni) contains t{N) 
and t{Ni) = t{N), e(iVi) = e(A^). N has an associated set of edges E{N) with {e{Ni) : 1 < i < 
k} C E{N) C -f{N). Let Ni denote the contraction ofV{Ni) in G. 

Case A^ ^ V: E{N) forms a (spanning) cycle on the vertices Ni, i = 1, ... , k. 

Case A^ = V.' Either (i) E{N) gives a cycle exactly as in the previous case, or (ii) E{N) is a 
spanning tree on the nodes Ni, rooted at Ni, with each e{Ni) the edge from Ni to its parent. 

Note that {e{Ni) : 1 < i < k} = E{N) in case (ii) above. In contrast when E{N) is a cycle, 
fewer than half its edges may belong to {e{Ni) : 1 < i < k}. See edge df in FigiJJ 

For any vertex v, a top-down traversal of T gives a naturally defined vt-path p{v) that starts 
with e{v), that we now describe. As an example in Figl2]p(j) = j, f,g,h,i,e,c,d,t; in FiglS] this 
path is composed of pieces in the subgraphs of 3 nodes, j, /; /, g, h, i, e; and e, c, d, t. 

For any interior node A'" let p{v, N) = p{v) j{V{N)). So p{v) = p{v, V). We will specify p{v) 
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Figure 3: Shortest-path structure. Node labels are z values. E{N) edges are the dashed edges 
joining the children of N. 

by describing the edge sets p{v, N). We leave it to the reader to add the simple details that specify 
the order of these edges in the desired vt-path. 

Consider any interior node N and a vertex v G V{N). The ends of p{v,N) are v and t{N). 
(This is clear for A^ = V, and we shall see it holds for the other nodes by induction.) Let v belong to 
V{Ni) for the child Ni of A^ (possibly v = Ni). E(N) contains a unique A^jA^i-path P that begins 
with the edge e{Ni). (For i = 1, P has no edges.) p{v, N) has the form 



piv,N) = E{P)u\J{p{. 



Xj , Nj J 



N, G V{P)}. 



(10) 



Implicit in (llOp is that p{v, Nj) has the form p{xj, Nj), i.e., it has t{Nj) as one of its ends. To prove 
this consider three cases: If Nj is the last node of P then p{v, Nj) ends at the end oi p(v, N), which 
is t{N) = t{Nj). If Nj is the first node of P then we have chosen the first edge of P as e{Nj), and 
it has t{Nj) as an end. If Nj is neither first nor last in P then P contains two edges incident to 
Nj, one of which is e{Nj). For what follows let fj be the edge not equal to e{Nj) that is incident 
to Nj. Observe that fj exists unless Nj is first in P. 

It remains to specify vertex Xj. This vertex is v if Nj is first in P, else it is the vertex fjriV{Nj). 

This completes the definition of p{v,N). Note that when N = V and E{N) is a tree, e{Nj) is 
the edge that leaves Nj in (the directed version of) a path p{v)] in all other cases (i.e., E[N) is a 
cycle) this needn't hold, e.g., e(/i) doesn't leave h in p{j). 

Next we specify the numeric values that will prove the p{v) 's are shortest paths. For any v &V lei 
Py denote a shortest vt-path. Our approach is based on the subgraphs PxUPyUxy. For motivation 
we first discuss our proof as it specializes to a nonnegative weight function. Recall Bellman's 
inequality for an edge yx, d{x) < d{y) + w{xy); rewrite it, with the above subgraph in mind, as 
2d{x) < d(x) + d(y) + w{xy). Consider an arbitrary vertex v. Assume equality holds in Bellman's 
inequality for each edge of yx \n P^, x ^ t. We show this implies w{Py) < w{P) for any t;t-path P. 
Let P = {v = xc^i, ... ,X£ = t), and add the inequalities d{xi) + d{xij^i) + i(j(a;jXj+i) > 2d{xi), to 
get 

d{v) + ^{2d{x) ■.xeP-v,t] + d{t) + w{P) > ^{2(i(x) : x G P - t}. 



Equality holds \i P = P^. Add in the identities 2d{x) = 2d{x) for each x ^ P to get 
d{v) + ^{2d{x) :xeV-v,t} + d{t) + w{P) > ^{2d{x) :xeV-t}. 
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Since equality holds for P = P^ we get w{P) > w{Pi,) as desired. The generalization to conservative 
weights will use a laminar family instead of the lower bounds 2d{x). 

A gsp-structure consists of a gsp-tree plus two functions d, z. Each vertex v has a value d{v), 
its distance to t. Each node N of T has a value z{N) that is nonpositive with the exception of 
z{V) which has arbitrary sign. Enlarge E{G) to the set Ei{G) by adding a loop xx at every vertex 
except t, with w{xx) = 0. Also for any such x define E{x) to be {xx} (although x is a node of T, 
its set E{x) has not been previously defined); set E{t) = 0. Say that a node N of T covers any 
edge with both ends in A'^ (including a loop xx at x £ N — t) as well as the edge e(A^) (if it exists). 
Every edge xy E Ei{G) satisfies 

d{x) + d{y) + w{xy) > y^{z{N) : N covers xy}, (11) 

with equality holding for every edge of |J{£^(A^) U e(A^) : N a node of T}- 

We show this structure guarantees that each path p{v) is a shortest vt-pat\i, by an argument 
similar to the nonnegative case above: Take any vt-path P~ . Enlarge it to a spanning subgraph 
P by adding the loops xx, x ^ V{P~). Similarly define Py to be a spanning subgraph formed 
by adding loops to p{v). Let cov{P,N) be the number of edges of P covered by N. Adding the 
inequalities (llip for each edge of P gives 



d{v) + '^{2d{x) ■.x£V-v,t}+ d{t) + w{P) > ^{cov(P, N)z{N) : A^ a node of T} (12) 

with equality holding for P^. 

Claim Py achieves the maximum value of cov{P,N), for every node N. Furthermore every P has 
the same value of cov{P,V). 

The claim implies the right-hand side of (112^ achieves its minimum when P = Py (recall z{N) < 
for every A^ 7^ V). This implies w{P) > w{Py) as desired. 

Proof of Claim: Every P and A^ satisfy cov(P, A^) < |A^ — t|. In proof, each vertex of A^ — t is 
either on a loop of P or on a maximal subpath of P from a vertex of P to an edge of 6{N) or to t. 
cov(P, A^) equals |A^ — t\ decreased by the number of the subpaths that leave A^ on an edge 7^ &{N). 
This also shows cov(P, A^) equals | A^ — t\ for P = Py and for A^ = V, as desired {> 

6.2 Basic facts 

This section derives the basic structure of shortest-paths to t for conservative graphs. The key 
concept, the "planted-cycle", is essentially a special case of the key concept for /-factors, the 2/- 
unifactors. Secj?] discusses how planted-cycles relate to general /-factor blossoms (see especially 
Lemma [28]) . 

For convenience perturb the edge weights by adding e* to the ith edge. Here e > is chosen 
small enough so that no two subsets of E^ have the same weight. Also assume the edges are ordered 
arbitrarily except that the loops xx are the last n edges. Let Py denote the unique shortest ft-path. 
Let Fy be its enlargement to a spanning subgraph, i.e., Py plus the loop xx for every vertex x ^ Py. 
(In the language of /-factors, Fy is the unique minimum /^-factor, see SecJT)) 

We start with a simple situation where conservative weights do not differ from general nonneg- 
ative weights. 
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Lemma 26. uv £ Pu — Pv implies Pu = u, v, P„. 

Proof. Let P^ = u,v,Qj], where Q^ is a vt-p&th, and assume Q^, ^ P^. This assumption imphes 
w{Qv) > w{Py), so the tit-trail T = u,v,Py has w{T) < w{Pu), whence T is not simple. The 
assumption uv ^ P^ implies T = C, Ru for a cycle C through u and a ut-path i?^. w(C) > implies 
w{Ru) < w{T) < w{Pu), contradiction. ■ 

This gives a weak analog of the shortest path tree: 

Corollary 27. {uv : uv £ Pu® Pv} is a forest. 

Proof. For contradiction let x^, x^, ... ,x^ be a cycle of these edges. Wlog assume x^x^ G P^-i — P.j.q. 
Lemma [26] shows P^-i = x^^^x^^P^o. This implies x^x'^ is not in P^-i so it is in Pr^2 — P^i. Thus 
Pr^2 = x^jX-'^jx'^, Pj-o. Continuing this way gives P^r = x^,x'^~^, ... ,x^,P^o. Px^ a path implies 
x^ ^ x^, contradiction. ■ 

We will consider subgraphs of the form P^ U P^, U uv, viewed as a multigraph contained in 2G. 
Define a p- cycle ("planted-cycle") to be the union of a cycle C and 2 copies of a path P from a 
vertex c £ C to t, with V{P) H V{C) = {c}. (Possibly c = t.) 

For motivation first suppose all weights are nonnegative. It is easy to see (e.g., using the 
shortest-path tree) that 

(a) uv ^ PuD Py implies P^ U Py U uv is a p-cycle. 

(b) uv £ PuC] Py never occurs (assuming our perturbation of w). 

In the remaining case uv £ Pu — Pv, the lemma shows P„ U P„ U uv consists of 2 copies of Pu- 
This can be viewed as a p-cycle if we add a loop uu to the graph. This motivates an approach 
similar to the algorithm of |14| for matching: Enlarge the graph by adding zero-weight loops xx, 
X £ V. Define the quantity C{uv) = w{Pu) + w{Pv) + w{uv). Repeatedly shrink a cycle of edges 
with minimum Q value (updating the weights of incident edges in a natural way). The cycles that 
get contracted are parts of the desired paths Pu. Once all these pieces are found we can assemble 
them into the complete paths. (The first edge in P^ is revealed in the step that shrinks the loop 

XX.) 

A more precise statement of this "shrinking procedure" is given in SeciT] (or see |14|): FigJH 
in SeclT] will show its execution on Figl2j Once stated, it is a simple exercise to check that the 
shrinking procedure gets the desired shortest path structure when weights are nonnegative. In fact 
the shrinking procedure becomes a variant of Dijkstra's algorithm. 

We wish to extend this to conservative weight functions. Observe that in Fig. [2] P^ U Pfe U a5 
is not a p-cycle - it contains an extra copy of ah. We can remedy this by deleting the extra copy. 
This suggests the following definitiono 

{w{Pu) + w{Pv)+w{uv) UV^PuHPv ,^^, 

Ciuv) = < (13) 

[w{Pu) + w{Pv) - W{UV) UV£PunPv. 

This definition has a similar failure: In Fig. [2] P^ U P/j U (^/i contains extra copies of gh, cd and ce. 
But we will see this failure is irrelevant to the shrinking procedure, and once again it constructs 
the desired shortest path structure. The reason is that C(5'^) = 3 is not the smallest (" value, and 



^These expressions are central for general /-factors - see (^uv and (^"" defined at the start of Sec l8.ll 
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FigJUwill show it gets "preempted" by C(ce) = 1 < 3. (Note also that property (a) fails in Fig. [2) 
Prf U P/ U df contains extra copies of cd and ce. Again FigJUwill show this is irrelevant.) 
For the rest of this section assume 

\m\ = 1- 

(This can always be achieved by adding a dummy edge incident to t. The shortest-path structure for 
the given graph can easily be derived from the structure for the enlarged graph.) The assumption 
implies any p-cycle has c^ t. 

Let cc' be the multiplicity 2 edge incident to C. We use the notation C, c, c' throughout the 
discussion. When convenient we do not distinguish between the p-cycle, or C, or the pairs C, c or 
C, cc' . Note that each vertex x £ C — c has exactly 2 xt-paths in the p-cycle, depending on which 
direction we traverse C. 

The following lemma gives properties of p-cycles that are shared with general blossoms, as well 
as properties that are specific to shortest paths (see SeclT]). 

Let E* be the set of edges in {uv : uv ^ P^® Pv} whose ("-value is smallest, and let Q* be this 
smallest (-value. 

Lemma 28. (i) There is a p-cycle whose cycle C consists of E* and possibly other edges of (^ -value 

< C- 

(a) For every x £ C, Px is one of the xt-paths in the p-cycle. In particular Pc starts with edge 



cc'. 



(Hi) Every x £ C has 2w{Px) < C* ■ 

Proof. First note the lemma is trivial if E* consists of a loop, so suppose not. The bulk of the 
argument consists of several claims that lay the foundation for {i)-{iii). 

Take any edge uv £ E* . Consider the subgraph S of 2G formed by modifying Pu U P^, to contain 
exactly one copy of uv. Note that 

w{S) = C{nv) = C. (14) 

We will show S is the desired p-cycle. 

Claim 1: E(S) can be partitioned into 

(a) a circuit C = {Pu © Py) + uv; 

(b) a multiplicity 2 path P from t to some c £ C (V{P) n V{C) = {c}); 

(c) zero or more multiplicity 2 paths joining two distinct vertices of C — c. 

Proof: We start by analyzing the multiplicity 1 edges of S. We work in the enlarged graph with 
edges Ei{G). Let Fu ^ Ei{G) be P„ enlarged to an /^-factor by adding loops; similarly for F^. The 
subgraph F^ ® F^ contains a nu-trail T that starts with an edge of F^ — F^ incident to u, ends with 
an edge of F^ — F^ incident to v, and has edges alternating between P„ — F^ and Fu — F^. (Any of 
these edges including the first and last may be loops.) P„ © T is an /j,-factor, so 

w{Fu) - w{Fu n T) + w{Fu n T) < w{Fu). 

Similarly 

w{Fu) - w{Fu n T) + w{Fu n T) < w{Fu). 

Adding these inequalities gives vu(Fu) -\- w{Fu) < w{Fi,) -\- w{Fu). Thus all inequalities hold with 
equality. The perturbed weight function implies Fu — Fu OT -\- F^, CiT = F^,. In particular Fu — T = 
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F^ — T. So the nonloop edges in T + uv are precisely the multiphcity 1 edges in S. Also we get (a) 
of the Claim. 

If T contains a loop xx, x ^ u,v, then the 2 edges that immediately precede and follow xx in 
T are both in P„ — P^ or both in P^ — P^- This implies Pu ® P,; is a trail consisting of 1 or more 
subpaths, each of 1 or more edges, that alternate between P„ — P^ and P^ — P^- 

S contains a multiplicity 2 path from t to some c ^T (with no other vertex of T) . We get (b) 
of the Claim. 

Let X be the set of all vertices of T — c that are on an edge oi P^C] P^ — uv. li x ^ u,v \s the 
end of an edge P„ n P^ — uv then x is on another edge of P^ and another edge of P^. It is easy to 
see this implies that each vertex of X is joined to another vertex of X by a multiplicity 2 path of 
S. These are the paths of (c) of the Claim. 

Finally note that (a)-(c) account for all edges of S: An edge of P„ © P^ is in C . An edge of 
Pu n P-u is in a path or cycle of such edges. A path either has both ends in C (making it type (c)) 
or one end in C (making the other end t, so the path is type (b)). A cycle of Pu H P^ edges cannot 
exist since Pu is acyclic. <^ 

The above X is a set, not a multiset, i.e., a vertex x G X is the end of only 1 type (c) path. 
This follows from d(x, S) = 4, which also holds if x = c. It is worthwhile to describe the case x = u: 
If uv & Pur\ Pu then d{u, S) = 2, so u ^ X. Suppose uv ^ Pu^ Pv and u £ X. The first edge of Pu, 
say /, must belong to Pu] let g be the other edge of Py H 5{u). Then S D 5{u) consists of 2 copies 
of / plus the edges g,uv. 

Claim 2: The edges of type (a) and (c) can he partitioned into a collection of cycles. 

Proof: Any circuit is the edge-disjoint union of cycles. So we can assume there are type (c) edges, 

i.e., the above set X is nonempty. 

\X\ is even, so the vertices of X divide the edges of C into an even number of segments Cj. 
Partition the edges of C into 2 sets Ci , C2 , each consisting of alternate segments Cj . Partition the 
edges of the type (c) paths into 2 sets Pi,P2) each consisting of 1 copy of each type (c) path Pj. 
Now the edges of type (a) and (c) are partitioned into the two sets nonempty Vs U C^, s = 1, 2. 

Let VC be one of the subgraphs Vs U Cg- Observe that each vertex x has d{x,VC) even: If x is 
interior to a Pj then d{x,VC) =2. If x is the end of a Pj then it is the end of a corresponding Cj, 
so again d{x,VC) = 2. (Recall that X is not a multiset!) Any other x has 1 or 2 occurrences in 
VC, both interior to Cj's. (2 occurrences may correspond to x occurring twice in some Ci, or once 
in two different Cj's.) Thus d{x,VC) E {2,4}. 

We conclude that each connected component of VC is a circuit (the construction ensures the 
circuit is edge-simple, since any edge of G occurs at most once in VC). So VC is a union of cycles. 

Claim 3: C is a cycle and there are no type (c) edges. 

Proof: Suppose the partition of Claim 2 consists of a single cycle. Then there are no type (c) edges 

(since each of the above sets Vg U Cs is nonempty). So Claim 3 holds. 

Now for the purpose of contradiction assume the partition of Claim 2 contains at least 2 cycles. 
Each cycle has nonnegative weight. In fact the perturbation implies the weight is positive. 

One of these cycles, call it B, contains vertex c. Form a subgraph S' by using the cycle B and 
the type (b) path of S from t to c. Our assumption implies w{S') < w{S). 

Let B contain an edge xy ^ Px®Py (Corollary [27]). Let B^ {B'^) be the xi-trail contained in S' 
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that avoids (contains) xy, respectively. Define By and B' similarly. Then 

w{B^) + w{By) + w{xy) = w{B'^) + w{B'y) - w{xy) = w{S') < wiS) = C- (15) 

(The last equation is (I14p .) The definition of Ci^u) shows it is at most either the first expression 
of (jlSp (ifxy ^ PjjUPj,) or the second expression (if xy G PxCiPy). Thus ^{xy) < C*, contradiction. <0 

Now we prove {i)-{iii). Take any edge xy E C. The relation (|15p becomes w{Bx) + w{By) + 
w{xy) = w{B'^) + w{B'y) - w{xy) = w{S) = C- So C{xy) < C- If xy^P^® Py we get C{xy) = C, 
and Px,Py is either Bx,By or B'^^B' i.e., [ii) holds for x. 

To prove (i) it remains only to show C contains every edge u'v' € E* . Analogous to ([H]), the 
subgraph S' formed from u'v' the same way S is formed from uv has weight w{S') = Q{u'v') = C,* . 
The perturbed edge weight function implies 5 = S' . Thus exactly 1 copy of u'v' belongs to 5, i.e., 
u'v' G C. 

To prove the first assertion of [ii) we need only treat the case xy ^ Px® Py. (Note the second 
assertion of (u)is a simple special case of the first.) These edges are a proper subset of C (Corollary 
[27j) . Let Q be a maximal length path of such edges that does not contain c internally. At least one 
end of Q, say r, is on an edge rs £ C — [Pr ® Ps) (the other end may be c). If rs ^ PrD Pg then 
Pr is the ri-path that avoids rs in C. Lemma [26] shows any x £ Q has Px a subpath of Pr- Thus 
{ii) holds for x. Similarly if rs £ Pr O Ps then Pr is the ri-path that contains rs in C, and Lemma 
[26] shows any x £ Q has Px the xi-subpath of C containing rs. Again (ii) holds for x. 

To prove (Hi) first assume x ^ c. Thus C* = w{S) = w{Bx) + w{B'^) > 2w{Px). For x = c the 
argument is similar: Since G is conservative, ^* = w{S) = 2w{Pc) + u^(C) > 2w{Pc). ■ 

We want the above cycle C to be a cycle node in the gsp-tree. Lemma [28t^ iz) shows C has the 
required properties for vertices x £ C. Now we show C has the required properties for x ^ C. 

A vertex x respects a p-cycle if either x £ C and Px is an xi-path in the p-cycle, or x ^ C and 
Px either contains no vertex of C or it contains exactly the same edges of 7(C) + cc' as some Py, 

Let C, cc' be the p-cycle of Lemma [28] 

Lemma 29. Every vertex x respects C . 

Remark: We allow C to be a loop cc £ Ei[G). In this case the lemma states that a shortest path 
Px that contains c actually contains cc' . 

Proof. Lemma [28l^u) shows we can assume x ^ C. The argument begins similar to Claim 1 of 
Lemma [28] We work in the graph with edges Ei{G). For any vertex u let F^ ^ Ei{G) be P^ 
enlarged to an /^-factor. The subgraph Ex ® Ec contains an xc-trail that starts with an edge of 
5{x) r\Ec — Ex, ends with an edge of 5{c) n Ex — Ec, and has edges alternating between Ec — Ex and 
Fx — Ec. Let e be the first edge of the trail that belongs to 6{C). Let T be the subtrail that starts 
at X and ends with edge e. 

Consider two cases: 

Case e 7^ cc' : Clearly e £ Ex — Ec. Let e be incident to vertex a £ C. T is a subgraph of 
Fx ® Fa (Lemma I28( u)). Now follow the argument of Lemma [28] Claim 1: Ex®T is an /^-factor 
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and Fa (BT is an /^-factor, so we get F^ — F^CiT + FaCiT = Fa and F^ — T = Fa — T. The latter 
implies Fx and Fa contain the same edges of 7(C) + cc' , i.e., x respects C. 

Case e = cc' : Clearly cc' € Fc — F^. Define a degree-constraint function f^ by 



fix) 



1 X = t 

3 X = c 

2 x^t,c. 



An /'^-factor (in Ei(G)) consists of a ct-path plus a cycle through c, plus loops at the remaining 
vertices. Since G is conservative and w has been perturbed, cc is the smallest cycle through c. So 
F'^, the minimum- weight /'^-factor, consists of P^ plus a loop at every vertex except t (in particular 
cc G F^). Thus F" = Fc + cc. 

Now the argument follows the previous case: T is a subgraph of F^. © F^. F^. ® T is an /'^-factor, 
F'^ e T is an f^-iactov, so F^ - F^ n T + Fc n T = F" and F^ - T = F" - T. We have already 
noted cc' ^ Fx and the last equation shows Fx contains the same edges of 7(C) as F^, i.e., every 
loop aa, a £ C. In other words Px does not contain a vertex of C, so x respects C. ■ 

6.3 Construction of the shortest-path structure 

The last two lemmas show how to construct the first node of the gsp-tree. We construct the 
remaining nodes by iterating the procedure. This sections first shows how to construct the gsp-tree; 
then it completes the gsp-structure by constructing z. 

For the gsp-tree we first state the algorithm and then prove its correctness. Shrinking a p- cycle 
means contracting its cycle C; C is the shrunken cycle. Let G be a graph formed by starting with 
G and repeatedly shrinking a p-cycle. (So the collection of shrunken cycles forms a laminar family.) 
We treat G as a multigraph that contains parallel edges but not loops. It is convenient to refer to 
vertices and edges of G by indicating the corresponding objects in G. So let E{G) denote the set of 
edges of G that correspond to (nonloop) edges in G. Thus writing xy G E{G) implies x,y G y{G). 
We do not distinguish between xy and its image in G. Similarly, writing C, cd for a shrunken p-cycle 
implies cc' G E{G), and a ft-path in G has v G V{G) and its first edge incident to v. An overline 
denotes quantities in G, e.g., w, (. 

The following algorithm constructs the gsp-tree T. We assume the shortest paths Px are known. 
(The function z is constructed below.) 

Initialize G (the current graph) to the graph (V,Ei{G)), and T to contain each vertex of G as 
a singleton subtree. Then repeat the following step until G is acyclic: 

Let C, cc' be the p-cycle of weight C*iG) given by Lemma[28l Shrink G in G. Set e(C) = cc' . 
Unless C is a loop, create a node in T whose children correspond to the vertices of C. 

When the loop halts create a root node of T whose children correspond to the vertices of the final 
graph G0 

To complete the description of this algorithm we must specify the weight function W for G. Let C 
be the collection of maximal shrunken cycles that formed G. An edge cc' G E{G) may be associated 



^"^In contrast with the general definition, our assumption that i is on a unique edge ensures the root of T is always 
a tree node. 
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with two cycles of C, one at each end. For a p-cycle C,cc' in C, as in Definition 1251 V(C) denotes 
the set of vertices of G that belong to C or a contracted vertex of C . ^(C) is the union of all the 
V{C) sets. For C,cGC and x € V{C), Cx ^ E[G) denotes the minimum- weight xc-path contained 
in7(C,G). 

Let G denote the given graph, and let 

u} = 2\w\{E). 

The weight of an edge e = xy £ E{G) in G is defined to be w{e) = w{e) + A(e, x) + A(e, y), where 



A(e,x) 



X ^ V{C) 

— 2lo C contains p-cycle C, e with x G y{G) 

2oj + w{Gx) C contains p-cycle C, cc' with x G V{C), cd ^ e. 



We will use a variant of the "respects" relation. Let S* be a cycle or path in G. Let C, cc' be 
a p-cycle of C, and C the contracted vertex for C in G. S respects C ii S H S{G) is either empty 
or consists of cc' plus < 1 other edge. This definition corresponds to the previous definition if we 
view C as a p-cycle whose cycle is a loop C G. More importantly we shall use this fact: If a path P 
respects G, cc' and G is derived from G by contracting C to a vertex G, then the image of P in G 
respects G. Also, for any C, cc' in C, if the image of a shortest path P„ respects G and contains G 
internally, then P^ traverses G along the path Gx that joins the 2 edges of Px n S{G). This follows 
from the optimality of P^- 

S respects C if it respects each G £ C. When S respects C, the preimage of S is the subgraph 
S of G that completes S with minimum weight, i.e., S consists of the edges of S plus, for each 
p-cycle C on 2 edges of S, the minimum-weight path Gx in G that joins the 2 edges. This preimage 
is unique. (Any C, cc' of C that is on 2 edges of S is on cc' and another edge that determines the 
vertex x in the definition of the preimage. Note that if 5 is a vt-path that respects C and v £ V{C) 
then V = c for some p-cycle G, cc' of C; so in this case the preimage of S does not contain any edges 
of C.) 

Lemma l32l below shows the following properties always hold. 

PI: G is conservative. 

P2: For any vertex x ^ ^(C), Px is the preimage of Px- For any p-cycle G,cc' £ C, Pc is the 
preimage of P^, both of which start with edge cc' . For any other x £ V{G), the image of Px in G 
is P^. 

Lemma 30. Assuming P1-P2 always hold, T is a valid gsp-tree for G. 

Proof. First observe that every vertex x £ V{G) gets assigned a value e(x). This holds as long 
as some iteration chooses xx as the minimum weight p-cycle. So for the purpose of contradiction, 
suppose an iteration contracts a p-cycle G where x £ G but the loop xx has not been contracted. 
Since xx ^ Px(BPx = 9, xx is considered in the definition of E* . So we get the desired contradiction 
by proving 

C(xx) = 'u;(xx) + 2w{Px) < C- (16) 

Lemma [28rizi) shows 2w{Px) < C* (iii ^W iteration). We can assume the perturbation of w gives 
every subgraph of 2G a distinct weight. Thus 2w{Px) < C* ■ Furthermore we can assume the 
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perturbation enforces a lexical ordering of the edges. Since the loop xx is ordered after any edge of 
G, it cannot reverse this inequality, i.e., (116p holds. 

The algorithm sets e(C) correctly for each node C of T, by definition. Recall the path p[v, N) 
from (fTOl) . It is easy to see the lemma amounts to proving that for any interior node C of T and any 
vertex x G V{C), the edges of E{Px) n 7(C) correspond to p{x,C) as defined in (fTO|) . Let p-cycle 
B, bb' be the child of C with x £ B. (B is a loop when the child of C is a leaf of T.) We consider 
two similar cases for node C. 

Suppose C, cc' is a cycle node of T- Let G be the graph immediately before C is contracted. 
P2 in G shows Pj, is the preimage of Pg-. It also shows Pg- starts with bb' . Lemma [28riz) shows Pg- 
is one of the B f-paths in the p-cycle. We have already observed that C and its children have the 
correct e- values. Thus 7~ specifies the desired path p{b,C) as defined in (|1U|) . This extends to any 
X G V{B) by the last assertion of P2. 

The remaining case is for the root node V of T, with V a tree node. In the final acyclic graph, 
P2 shows Pfe is the preimage of P-g-. Thus T specifies the desired path p{b, V). As before this extends 
to any vertex x £ V{B). (} I 

We complete the gsp-structure by specifying z. For each node N of T let ^jy be the weight 
C*iG) of its corresponding p-cycle. For a leaf x the corresponding p-cycle is the loop xx; for the 
root V of T, which is a tree node, the corresponding p-cycle is the last p-cycle to get shrunk. Let p 
be the parent function in T- For each node A^ define 

Lemma [32] shows this additional property always holds: 
P3: For any uv eE{G), 



C{uv) 



'({UV) UV<^Pu®Pv 

C{uv) uvePu- Pv,U^V{C) 

^C{uv)-Aio UV e Pu- Pv,u^V{C). 



The notation in P3 is unambiguous since uv G P„ iff uv G Pu by P2. 

Lemma 31. Assuming P1-P3 always hold, T with the above function z is a gsp-structure for G. 

Proof. The definition of z clearly implies that for any node N of 7~, the z- values of all the ancestors 
of N (including N) sum to C,n- We claim C*{G) increases every iteration. The claim implies that 
for every A^ / V, 2:(A^) = C,N — Cp{N) ^ 0> i-^-) •^(A^) is nonpositive as desired. To prove the claim 
note that for a fixed edge uv, uv belongs to P„ © P^ in one iteration iff it does in the next iteration, 
as long as it is not contracted (by P2). Hence the only change in E* from one iteration to the next 
is that contracted edges leave E* . P3 shows the edges in E* retain their original (^-values. Thus 
C*iG) never decreases. 

To complete the proof we must show (jlip for every uv G Ee{G), with equality for e(A^) edges 
and other edges of subgraphs E{N). Consider two cases. 
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Case uv ^ PuH P^ : By definition 

d{u) + d{v) + w{uv) = C{uv). 
Let A^ be the deepest node of T that covers uv. Observe that 

C{uv) > Cn, 

with equality holding if uv G E{N) U e{N). In proof, if uv ^ P^^ Py this follows from uv £ E*. If 
uv £ Pu — Pv, then Lemma [26] shows uv = e{u). So the loop uu covers uv, N = {u}, and equality 
holds. Since the sum of the right-hand side of (|lip equals (]\r, combining the two displayed relations 
gives the desired conclusion for (jlip . 
Case uv £ PuCi Py : By definition 

d{u) + d{v) — w{uv) = C{uv). 

uv is covered by the nodes of T that are ancestors of u or v. The ancestors of u have z-values 
summing to C(''^''^) = 2d{u) and similarly for v. Let node A be the least common ancestor of u and 
V in T- Since uv is in the cycle of node A, (^j^ = C('U'v) = d{u) + d{v) — w{uv) (by definition of Q. 
The sum of the right-hand side of (jlip equals 

Cu + Cv- C.A = '^d{u) + 2d{v) - {d{u) + d{v) - w{uv)) = d{u) + d{v) + w{uv) 

as desired. ■ 

The development is completed by establishing P1-P3: 

Lemma 32. P1-P3 hold in every iteration. 

Proof. Consider two edges Cj = Xiyi G E{G), i = 1,2, where both Xj have the same image in G. If 
Xi^V{C) then 

A(ei,xi) + A(e2,X2) = 0. (17) 

Suppose the Xi belong to the p-cycle C of C. If ei, 62 respects C then wlog the p-cycle corresponds 
to C, 62, and 

A(ei,xi) + A(e2,X2) = w;(C^J. (18) 

In the remaining case, i.e., 61,62 does not respect C, 

A(6i, xi) + A(e2, X2) > 2(w + \w\{E - 7(C))). (19) 

Claim 1: Let S he a cycle of G. If S respects C then w{S) = w{S) for S the preiniage of S. If S 
does not respect C then w{S) > 2uj. 

Proof: Let CTZ (CM) contain the p-cycles of C that are respected (not respected) by S, respectively. 
Let CX = U{Cx : x = xi in I^Sii for G G CTZ}. Then ([I7])-([l9]) imply 

w(S) > w{{e : 6 G 'SuGX}) + J^{2(w + \w\{E - j{G))) : G G CAf}. (20) 

Here we use the fact that the 2a; terms cancel on the 2 edges incident to a cycle of CTZ (this property 
actually depends on our assumption that t has a unique incident edge). 
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When S respects C, i.e., CM = 0, (|20p holds with equahty, and we get the claim. Suppose 
S does not respect C. Then choosing B as any p-cycle of CM, (|20p implies w{S) > w{{e : e E 
SUCX}) + 2{u + \w\{E - jiB))) > 2uj, giving the claim. <^ 

Claim 1 implies property PI. 

Claim 2: Let S be a vt-path in G. If S does not respect C then w{S) > lo. If S respects C let S be 
the preimage of S. If v ^ ^(C) then w{S) = w{S). If v £ V{C) then w{S) = w{S) — 2uj. 
Proof: S satisfies a version of (f20|) that accounts for the term A(ei,'y) for edge ei = 5{v,S). We 
examine several cases. 

Suppose S respects C. If w ^ ^(C) then (|20l) holds with equality, giving the claim. If -y G ^(C) 
then C contains a p-cycle C, ei with v G V{C). w{S) contains an extra term A(ei,f) = —2uj, again 
giving the claim. 

Suppose S does not respect C. If f ^ ^(C) then (|2Up holds unmodified. As in Claim 1, 
w{S) > 2a;, giving the current claim. So suppose v S V{B) for B £ C. If ei does not respect B 
w{S) contains an extra term A(ei,f) = 2uj + w{By) > uj + \w\{E — ^{B^)), giving the claim. In 
the remaining case ei respects B and S does not respect some A ^ B. The right-hand side of (|20p 
contains the extra term A(ei, t;) = —2a;, and the term for A is at least 3a; + 2|u;|(£' — 7(A)). These 
two contributions sum to > uj+\w\{E — j{A)). Thus the right-hand side of (j20p is > u, as desired. {> 

Claim 3: P2 holds every iteration. 

We argue by induction on the number of iterations. Consider an iteration for the p-cycle C, cc' . 
Let G~ ,w~ [G'^jW'^) be the graph and weight function immediately before (after) C is contracted, 
respectively. For greater precision, if H is G^ oy G^ and z is a vertex of H, let P{z, H) denote the 
shortest zt-path in H . Take any x G V{G). Let P = P(x, G^). Note that P^ continues to denote 
the shortest xt-path in G, and for P2 we want to establish the relationship between P and Px- C 
denotes the family of contracted vertices in G"*", i.e., it includes C . 

First assume either x ^ V{C) or x = b for some p-cycle B, bb' of C (x = c is a possibility). Let 
(5 be (2a;) if x ^ ^(C) (x = b) respectively. Let Q = P{x,G~). Lemma [29] shows Q respects C. 
Thus Q has an image Q'^ in G "*" that respects G. The inductive assumption of P2 in G ~ shows 
Px is the preimage of Q. (For x = c, x in G ~ may be a vertex or a contracted cycle, and we use 
the appropriate assertion of P2.) Thus the optimality of Px implies it is the preimage of Q^. Thus 
Claim 2 shows 

w~{Q) = w+{Q+)=w{Px)-6. 

Since a; > w{Px), Claim 2 shows P respects C and 

w+{P) = w{p-)-5 
for P~ the preimage of P. Since P~ is an xi-path, 

W{P-) > w{Px). 

Combining the inequalities gives w^{Q^) < w^{P). Thus Q^ = P. So as asserted by the first 
part of P2, Px is the preimage of Px- 

Also, as in the second assertion of P2, every p-cycle B, bb' of C has Pf, and P{B, G~^) = P{b, G ~^) 
both starting with edge bb'. This follows since P(b, G ~) starts with bb' (by definition if b is not in a 
contracted vertex of G ~, else by P2 in G~) and Pj, is the preimage of both P{b, G~) and P(6, G "*"). 
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Finally the last assertion of P2 follows since any y G y{B) has y = b. 

Claim 4: P3 holds in every iteration. 
Proof: By definition 

((uv) = w{Pu) + w{Pv) e w{uv) = {w{Pu) e A{uv, u)) + {w{Pv) e A{uv, v)) e w{uv) 

for some © G {+, — }. Since uv is an edge of G, P2 shows uw E P^ iff ^^ G Pu- Hence the version 
of this formula for Cl"*^"^) uses the same sign ® as C(^^)- We evaluate ({uv) using the formula, with 
Claim 2 providing W{Pu) and the definition of A{uv,u) giving its value, as follows. 
Suppose u ^ ViC). Pu is the preimage of Pu by P2. So 

wiPu)®A{uv,u) = w{Pu) 

since W{Pu) = w{Pu), A{uv, u) = 0. 

Suppose u E ^(C), say u E y{C) for the p-cycle C,cc'. Pc is the preimage of Pu by P2. Claim 
2 shows w{Pu) = w{Pc) — 2uj. If uv ^ Pu then 

since A{uv, u) = 2u} + w{Cu), and P^j = C^ U Pc by P2 and the optimality of Pu- If uv E Pu then 

w{Pu) - A{uv,u) = Uj{Pu), 
w{Pu) + A(n7;, u) = w{Pu) - 4w 

since P2 implies u = c and A{uv,u) = —2uj. 

The alternatives of P3 all follow, by combining the equations with — signs for uv £ PuD Pu and 
+ signs in all other cases. () I 



7 Background on /-factors 

The bulk of this section reviews the approach of |14| based on critical graphs. The review ends by il- 
lustrating the shrinking procedure. Then we show that procedure efficiently constructs a generalized 
shortest-path structure. 

To define critical graphs, for each vertex v £ V define /„, the lower perturbation of / at v, by 
decreasing f[v) by 1. Similarly define /", the upper perturbation of / at f, by increasing f{v) by 1. 
Every /„,/", t; E ^ is a perturbation of /. f'l.v stands for a fixed perturbation that is either /„ or 
P . A graph is f -critical if it has an /'-factor for every perturbation /' of fPH 

It is easy to see that a maximum /-factor of G can be found by working on the critical graph 
G^ formed by adding a vertex s with edges su , v E y{G) and loop vv, all of weight 0, and /(s) = 1: 
A maximum /^-factor is the desired /-factor. 

Recall the linear programming formulation for maximum weight /-factors, especially the dual 
problem (e.g. |27l Ch. 32]). We summarize the slight modification used in |14| for critical graphs. 
Simply put we want the dual variables to be optimum for every perturbation of /. 



^^The usual notion of criticality for matching only assumes existence of the factors for lower perturbations. It is 
easy to see this implies existence for all the upper perturbations. This holds for &-matching too, but not for general 
/-factors [Ti] . 
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Dual variables for a graph with integral weights are functions y : y — )• Z and z : 2^ x 2^ — )■ Z. 
Pairs B = (y{B),I{B)) with nonero z-value are called (weighted) "blossoms". z{V) is a shorthand 
for 2;((y, 0)). The optimum dual function y is defined by y{v) = —w{Fy) |14| Theorem 4.17], as 
expected from SeclHl 

The z function has support given by a forest that we now describe. It generalizes the shortest- 
path structure of Section [HI The usual version corresponds exactly to that structure; it also cor- 
responds to the blossom tree for matching |14| . Our algorithms use a weighted version of this 
structure - it is the shortest-path structure/blossom tree with every cycle node/blossom of z 
value contracted into its parent. This weighted blossom forest is defined as follows (we give a minor 
modification of |14| Definition 4.14], using |14| Lemma |141 Lemma 4.12]]). Let Ai be the set of 
all maximal blossoms of G. M = {V} for matching and 6-matching but not generally (for shortest 
paths M is the set of maximal cycle nodes) . 

(i) Each B £ Ai is the root of a weighted blossom tree W. Each interior node of W is a weighted 
blossom and each leaf is a vertex of G. The children of any node B are the maximal weighted 
blossoms properly contained in B plus all vertices of G contained in B but no smaller weighted 
blossom. V{B) is the set of all leaf descendants of i?, i.e., the vertices in B. Each vertex of G 
belongs to exactly one blossom of Ai . 

[ii) The support of z is {B, V : B a blossom of a W-tree}. z{B) > for blossoms B of W, 
while z{V) may have arbitrary sign. 

(iii) Each blossom B has a set of edges I{B) C 5{V{B)). 

[iv) The blossoms of Ai are nodes of a tree T. Every edge AB of T belongs to I {A) © I{B)- 

For shortest paths, I{B) consists of the edge e{N). As in shortest paths and matching, a blossom 
B covers any edge of ^[V{B)) U I{B). An edge uv of G has a value 

yzie) = y{e) + z{B, V : B a. blossom of W that covers e}. 

(Here we use the above convention; recall e = {u,v}.) e is underrated if 

w{e) > yz{e); 

e is strictly underrated if the inequality is strict and tight if equality holds. An /;|:i'-factor F respects 
a blossom B iff 

Fn<5(y(5)) = | ^Sf] ""^I^Sm ■ J • u^rfT,^^ (21) 

"- "^ " y 1(B) © e V ^ V{B), e is some edge m 6{V{B)) ^ ' 

(For shortest paths, the first alternative says a shortest ui-path leaves B on e(i?) if v € V{Ji)\ 
otherwise the second line says it contains e(i?) and one other edge (t ^ B) or one edge ^ e(^) 
(t G -B)). An /;|:v- factor has maximum weight if it contains every strictly underrated edge, its other 
edges are tight, and it respects every blossom with positive z-value. (The optimum /-factors we use 
satisfy this criterion.) 

As in matching, blossoms are built up from odd cycles, defined as follows [141 Definition 4.3]: 
An elementary blossom B is a 4-tuple {VB,G{B),GH{B),I{B)), where VB C V, C{B) is an 
odd circuit on VB, CH{B) C -i{VB) - E{C{B)), I{B) C 5{VB), and every v eVB has f{v) = 
d{v,G{B))/2 + d{v,CH{B)Ul{B)). G{B),GH{B), and I{B) the circuit, chords, and incident edges 
of B, respectively. We sometimes use "blossom" or "elementary blossom" to reference the blossom's 
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Figure 4: Contracted graphs for nonloop unifactors. 



circuit or its odd pair. Blossoms for shortest pairs have |/(5)| < 1 (as mentioned) and CH[B) = 
(this is the import of Lemma [28l especially part [ii) and Claim 2 respectively). 

The blossoms for optimum duals are found by repeatedly finding the next blossom and shrinking 
it. Fig m illustrates the shrinking of the 3 cycle nodes of FigJSl e.g., {a, 6, c} gets shrunk in Figj^a). 
In general the vertices of the blossom are contracted and a loop is added to the contracted vertex. 
Edges incident to the contracted vertex get their weights adjusted to account for contracted edges. 
(For shortest paths these edges are the paths p{v,N). In FiglH^a) w{be) decreases by 1 to account 
for path b, a, c.jHI 

For simplicity perturb the edge weights slightly so that no two sets of edges have the same 
weight. That is, number the edges from 1 to tti and increase the weight of the ith edge by e* for 
some e > 0. For small enough e > 0, no two sets of edges have the same weight. Thus any such 
perturbation has a unique maximum factor which is also maximum for the original weights. Let F^ 
and F'" denote the maximum weight fy and f" -factors respectively. (Eventually (in ()22p ) we set e 
to and define the dual function z using the original weights.) Recall (Sec. [2]) the multiset notation 
2S, S/2, 2G. Assume that even in 2G, no two sets of edges have the same weight. 

We choose the blossom to shrink next using a subgraph that generalizes the u7;-subgraphs of 
Section (6) Any subgraph of 2G consists of edges of G at multiplicity 0,1 or 2. A 2f-unifactor is 
a 2/-factor of 2G whose multiplicity 1 edges form an odd circuit. (When possible we abbreviate 
"2/-unifactor" to "unifactor".) The elementary blossom for a unifactor U with odd circuit C is the 4- 
tupleS = {V{G),G,CH{B),I{B)) whei'eCH{B) = {Unj{C)- E{G))/ 2 and I {B) = {Un6{G))/2. 
For shortest paths note how this defines I{B) to be e{N) (the type (b) path to b ends in multiplicity 
2 edge e{N)). 

The shrinking procedure of |14| constructs the blossoms as follows: 



Let the next blossom B be the elementary blossom of the maximum weight proper 2/- 
unifactorP^I Shrink B. Repeat this step until no unifactor exists. 



^■^The weight adjustment of [141 Fig. 13] also involves a large quantity J. We omit it since J is needed only for the 
proof of [14| . not for explaining the result. 

'^■^"Proper" means (a) the unifactor's circuit is not a loop created previously when a blossom was contracted; (b) the 
unifactor respects each previous blossom. Here "respects" is the generalization of (PT|) to unifactors, e.g., a shortest 
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Let us describe how this procedure constructs the blossoms and tree of Figl3j We use FigUl 
The shrinking procedure repeatedly finds the minimum weight proper unifactor. Let (^* denote its 
weight. C* never decreases from step to step, so we describe the blossoms found at the various values 
of C*- Recall that we treat the shortest-path problem by adding a weight loop at every vertex 7^ t. 

C* = —2: The loops at a, b and j are blossoms. (For instance the unifactor for j corresponds 

to the shortest j't-path - it consists of loop jj (multiplicity 1) and the jt-path of weight —1 

plus a loop at every vertex not on the path (all multiplicity 2). The weight of the unifactor, 

—2, is drawn inside the loop at j in FigHJ^a). Figl4]does this for all blossoms.) 

C* = 0: The loops at d, e, g, t are blossoms, as is cycle a, b, c. FigjJl^a) shows the graph after 

all these blossoms have been shrunk. 

C* = 1: Contracted vertex {a,b,c} plus vertices d,e,t form a blossom. FigHJ^b) shows the 

graph after it has been shrunk. 

C* = 2: The loops at /, h, and i are blossoms. 

C* = 3: All vertices but j form a blossom. FiglDj^c) shows the graph after it is shrunk. The 

new weight on jt reflects the contracted edge fd. 

Note that when a loop xx becomes a blossom, the unifactor 's weight is 2d{x). This corresponds 
to the bound 2d{x) in Secl6]mentione for the Bellman inequality argument. Also, the description for 
loops in C* = — 2 might seem to imply we need to know the shortest paths to execute this procedure 
- but see the implementation in SeclH] 

The optimum dual function z is defined as follows. For any blossom B let U{B) be its corre- 
sponding unifactor. If B is not a maximal blossom let p{B) denote the blossom B gets contracted 
into. Let w be the original (unperturbed) weight function. 

, , J w{U{B)) S a maximal blossom, , , 

^^^ " \ w{U{B)) - w{U{p{B))) B nonmaximal. ^ ' 

For instance in Figj3]blossom aa is created at C,* = —2, {a, 6, c} is created at C,* = 0, and z({a, b, c}) = 
-2-0 = -2. 

The following characterization of maximum weight unifactors is central to the analysis: 

Lemma 33 ( |14| Lemma 4.5 and Cor. 4.6]). For any vertex v, the maximum weight 2f -unifactor 
containing v in its circuit is Fy + F^ . For any edge uv consider 3 cases: 

uv ^ FuD Fi,: The maximum weight 2 f -unifactor containing uv in its circuit is Fu -\- F^ -\- uv. 

uv £ Fuf] Fy-. The maximum weight 2 f -unifactor containing uv in its circuit is F^ + F'" — uv. 

uv £ Fy — F^: Fu = F^ — uv. Furthermore the maximum weight 2f -unifactor containing v in 
its circuit is Fy -\- Fy -\- uv, and this unifactor contains uv as an edge incident to its circuit. 

For shortest paths (i) - (ii) characterize the maximum nu-subgraph. (Hi) corresponds to Lemma 



We close this review by reiterating some notation from |14) that we use in the next two sections: 
w: the given, unperturbed, weight function. 

B: the forest whose nodes are the elementary blossoms found by the shrinking procedure, plus 
the vertices of G. The children of a blossom B are the vertices that get contracted to form B. 



path that enters and leaves a blossom must contain e{N). [14] enforces (b) using the previously mentioned quantity 
J that we omit. 
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z{B): the sum of the dual values z{A) for every ancestor A oi B in B. If B does not appear 
in the weighted blossom forest W (defined above) then z{B) = z{A) for the blossom ^ of W that 
absorbs B by contraction. 

By, for any vertex v G V{G): the smallest weighted blossom containing v. 

8 Weighted /-factors 

This section presents the reduction of maximum weight /-factors to unweighted /-factors. As usual 
we use the shrinking procedure to find the weighted blossoms B in order of decreasing value z{B). A 
blossom is a pair of sets {y{B), 1(B)). We find these pairs in two steps using a graph (?(("): Having 
found the blossoms of z- value > C^, we construct G((). Its 2-edge-connected components constitute 
the T/-sets for all weighted blossoms of z- value (■ The edges of /-sets that are still unknown are 
found amongst the bridges of G(C) or in a related computation. 

In our references to |14| we are careful to recall that |14| modifies the given weight function w 
in two ways: First, the given edge weights are perturbed so every maximum factor F^v and every 
maximum 2/-unifactor is unique. Second, in |14| each time a blossom is contracted the weights 
of its incident edges are modified (by adding large quantities, including a value called W greater 
than the sum of all previous edge weights). The reduction of this section has no access to these 
conceptual modifications, and it must works entirely with the given edge weight function w. 

Our overall strategy is similar to the 6-matching algorithm of Section [5l We assemble the desired 
maximum weight /;|:D-factor from its subgraphs that lie in the various nodes of the weighted blossom 
tree. The details are similar to Section [5l and it is implicit in |14| that these details work. But for 
completeness the next lemma proves the necessary properties. 

Any blossom S in a W-tree has an associated graph G{B). Its vertices are the children of 
B, with every blossom child contracted. Its edges are the underrated edges that join any two of 
its vertices (this includes underrated loops). Informally the lemma states that the edges of any 
maximum /;|:i)-perturbation in G{B) satisfy all relevant constraints. 

Lemma 34. For any v £ V{G), any perturbation fl^v, and any B in W, let B C 5{B) be a set of 
edges that respects B. If B has the form / © e assume e is tight. Then G{B) U B has a subgraph F 
containing B wherein 

(a) every child of B that is a vertex x £ y{G) has fX^ix) = d{x,F); 

(b) F respects every child of B that is a contracted blossom. 

Proof. The two possible forms for B are 1(B) and 1(B) © e. We first show that 1(B) is unique and 
the lemma holds when B = 1(B). 

For any v G V(B), any maximum /;|:f -factor F respects B, i.e., it has F n 5(B) = 1(B). (This 
follows from the optimality of the duals.) So 1(B) is unique. Also this factor F proves the lemma 
when B = 1(B), i.e., F satisfies (a) and (b). 

For V ^ y(B) any set respecting B has the form 1(B) © e, e € S(B). By assumption e is 
tight. Let e = xy with y G V(B). If e G 1(B) then F = F^ — xy is a, maximum /a;-perturbation, 
with F n d(B) = I — e. (This follows from the optimality criterion of the duals.) Similarly if 
e G S(B) — 1(B) then F = Fy + xy is a, maximum /^-perturbation with F n 5(B) = I + e. As before, 
the maximum factor F proves the lemma. ■ 
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8.1 Finding the V{B)-sets 

The reduction is given the quantities, w{Fy),w{F^), v G V{G). So it can use these quantities: 

For each vertex v £ V, Cv = w^F^) + w{F^). 

For each edge uv G E, (uv = w{Fu) + w{Fi,) + w{uv) and C™ = 'w{F^) + w{F^) — w{uv). 

For a vertex v let UNI{v) be the maximum 2/-unifactor containing w in its circuit. Lemma [33] 
shows C„ = w{UNI{y)). Recalhng the definition of dual variables ([22]) we get 

Cv = z{By). 

For any edge uv of G let UNI{uv) be the maximum 2/-unifactor containing uf in its circuit. 
We classify uv as type 0, 1 or 2 depending on the number of sets F^, F^ that contain uv. Specifically 
uv £ E is 

type if uv ^ F^U Fy] 
type 2\iuv £ FuDFy-, 
type 1 \i uv £ Fu® Fy. Additionally a type 1 edge uv is type 1^, if C«) Cv > C- 

Of course the type of an edge is unknown to the reduction! Lemma [33] extends to give the following 
combinatoric interpretations of (uv and (^^'". In contrast to (^y these interpretations are also unknown 
to the reduction. 

Lemma 35. Consider any edge uv of G. 

(i) If uv is type then UNI{uv) = F^ + Fy + uv. Thus w{UNI{uv)) = Quv Furthermore 

(a) If UV is type 2 then UNI{uv) = F^ + F^ — uv. Thus w{UNI{uv)) = C"^. Furthermore 
/■uv ^ /- /■ <^ /- 

S -^ SUtSV -^ SUV 

(Hi) Suppose uv is type 1 with uv £ Fy — F^. Then UNI{v) = Fu + Fy + uv and uv is incident 
to the circuit ofUNI{v). So w{UNI{v)) = Cv = Cuv Furthermore C,^'" = C,u- 

Proof. Lemma [33] gives the characterization of the various unifactors. The relations between the 
various C, quantities all follow easily from this simple identity: Any edge uv £ E satisfies 

C™ + C™ = a + C.- (23) 

For instance to prove part (i), the relation w{Uuv) ^ 'w{Uu),w{Uy) translates to (^^v < CuXv- Now 
([23]) implies C«,C«<C''- ■ 



We define the graph G(C), for any real value C: Its vertices are the vertices of G with all blossoms 
of z-value > C contracted. Its edge set is 

E{GiC)) = {uv : uv G E{G), mm{Cuv, C™} > C}- 
The reduction can construct G{Q, since previous iterations have identified the blossoms of z- value 

Recall a blossom i? of fi is constructed as an elementary blossom in a contracted graph G. As 
such it has a circuit G{B). G{B) consists of edges that are images of edges of G, as well as blossom 
loops (resulting from contractions). In the lemma statement below, E{G{B)) n E{G) denotes the 
edges of G whose images belong to G{B). 
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Lemma 36. For any blossom B ofB, E{C{B)) n E{G) C E{G{z{B))). 

Proof. Consider any edge uv E C{B). The desired relation uv G E{G{z{B))) amounts to the 
inequality min{<^„„, ^"*'} > z[B). (Note that uv cannot be a loop in G{z{B)) since uv joins distinct 
blossoms of G{B).) We will prove the desired inequality using the fact that every edge uv € G{B) 
is tight, proved in |14| Theorem 4.17]o Consider three cases: 

Case uv is type 0: Tightness means y{u) + y{v) + z{B) = w{uv). Equivalently z{B) = w{Fu) + 
w{F^) + w{uv). Thus z{B) = C^^ and Lemma \35\i) shows uv G G{z{B)). 

Case uv is type 2: Tightness means y{u) + y{v) + z[Bu) + z{By) — z{B) = w{uv). As noted 
above 

w{F,) + w{F'') = Cv = z{B,) 

and similarly for u. Sustituting this relation gives w{F'^)+w{F^') — w{uv) = z{B). Thus z{B) = C^^^ 
and Lemma [351 m) shows uv G G{z{B)). 

Case uv is type 1: Wlog uv G F^ — F^. Tightness means y{u) + y{v) + z{B^) = w{uv). 
Equivalently 

z{B„) = w{Fu) + w{F^) + w{uv). 

With z{Bi,) > z{B) (since B is an ancestor of B^,) this gives ^„„ > z[B). 

The (last) displayed equation is equivalent to w{Fv) + w{F'") = w{Fu) + w{F^) + w{uv). Re- 
arranging gives w{F^) + w{F'") — w{uv) = w{F'^) + w{Fu) = z{Bu) > z{B), i.e., C™ > z{B). 



In G{C) every type edge is in a blossom circuit (Lemma ISST i). which refers to the given 
graph, and |14| Lemma 4.13], which shows the relation of Lemma [35li) is preserved as blossoms are 
contracted). The same holds for every type 2 edge. (However note that an arbitrary type or 2 
edge needn't belong to a blossom circuit - it may not appear in any G{C) graph because of blossom 
contractions.) 

In contrast a type 1 edge of G(C) may or may not be in a blossom circuit. The 1^ edges obey 
the following generalization of Corollary [27] for shortest paths. 

Lemma 37. In any graph G{Q the 1(^ edges are acyclic. 

Proof. We start with a relation between the set I{A) of a blossom A & B and the same set when A 
is contracted, i.e., set I (a) for blossom vertex a. In G{Q suppose a is a blossom vertex and vertex 
b y^ a. Recall that Fa Pi 5{a) = I (a) and Fi, respects blossom a ([Ml Corollary 4.11], which says 
that as expected, every maximum perturbation F'l^v respects every maximum proper unifactor's 
blossom). Thus 

\{Fa n 6{a)) e (Ffe n 6{a))\ = \I{a) (F^ n 6{a))\ = 1. (24) 



We use (I24p to prove the following: 

Claim In G(C) consider a blossom vertex c and distinct edges e, / G S{c). Suppose vertices b,c,c' 
have e£ Fb® Fc and f G FcB Fc'. Then f G Fh®Fc'. 

Proof: We can assume b ^ c since otherwise the lemma is tautologous. Since Ff, respects blossom 
c and e, / G 5(c), e G F^ (B Fc implies / ^ Fh ® Fc (by (p^ ). Combining with the hypothesis 



^*In matching and &- matching tightness is forced by the fact that every edge of a blossom circuit is in some maximum 
perturbation. The analogous statement fails for /-factors - an "exceptional" circuit edge may belong to no maximum 
perturbation at all or to every maximum perturbation [14| . 
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f (^F^(B F,, gives / G (Fb © F^) ® {F^ ® F,,) = Fb®F,,. 

Now consider a cycle of 1^^ edges in G{Q), say a, b, ... ,y,z with z = a. Type 1 means ab G Fa®Fb. 
This immediately shows the cycle has > 2 edges. We will show edge yz = ya also belongs to Fa®Fi,. 
ya may be parallel to ab (a length 2 cycle) but ya is not the same edge as ab (i.e., we do not have 
a length 1 cycle, since a loop aa is not type 1). So we get 2 distinct edges in 6{a) D {Fa ® Fh). This 
contradicts (|24p . 

It is convenient to also denote the cycle as a,c^ = b,c^, ... ,d'~^ = y,c^' = z. Type 1 means 
d~^c'^ E F^t-i © F^i. Inductively assume d~^d E F^ © F^i. (This holds for i = 1.) Since dd^^ E 
F^i © -Fc'+i the claim (with c = d, c' = d'^^) shows dd^^ E -Ffe © F^i+i. Thus induction shows 
d~^d E Ffe © Fcr, i.e., yz G Fb® F^, as desired. ■ 

The reduction processes the graphs G{C) for <^ taking on the distinct values in 

= 1 min{(^„^, (^"''} : uv an edge of G} 

in descending order. For any C E fi let (^~ be any value strictly between (^ and the next largest value 
in ri. Observe that G(C^) is the graph G(C) with every blossom of z- value > C contracted. 

Corollary 38. For any C ^ ^, G{C~) is a forest. 

Proof. In G(C) every type or 2 edge is in a blossom circuit (as indicated after Lemma [36l) . So 
G(C~), which has all these blossoms contracted, has only type 1 edges uv. Furthermore since 
nim{(uv, C"''} > (^, Lemma \35\ iii) shows Cu,Cv ^ C- ^^ other words uv is type 1a-. Thus Lemma 
[371 shows G(("~) is acyclic. ■ 

The next lemma shows how the reduction finds the vertex sets of the blossoms of W that have 
z- value (^. 

Lemma 39. As vertex sets, the 2- edge- connected components of G{C) are precisely the weighted 
blossoms of z -value (". 

Proof. We first show that each blossom of B is 2-edge-connected. More precisely let i? be a node 
of B with z{B) = C. 

Claim In G(C) the subgraph of edges [j{C{A) : node A of B descends from B and z{A) = (^} is 
2- edge- connected. 

Proof: The argument is by induction on the number of descendants A. Recall B is constructed by 
repeatedly finding the next elementary blossom B and contracting it. 

When B is found, each vertex of its circuit is either (i) an original vertex of G, or [ii) a con- 
tracted blossom of z-value > C,, or {Hi) a contracted blossom of z-value C,. (Recall the definition 
of z, (|22p . The type [Hi) blossoms are blossoms with 2- value 0.) Vertices of type (i) or [ii) are 
vertices of G(C)- Vertices of type [Hi) have 2-edge-connected subgraphs in G{C,) by induction. Each 
original edge of G{B) is contained in G{C,) (Lemma [36]). Since G{B) is a circuit when B is formed, 
it completes a 2-edge-connected subgraph of G{C,). This completes the induction. <^ 

Now starting with G{C,), contract each of the above 2-edge-connected subgraphs that corresponds 
to a maximal blossom of z-value Q. We get the acyclic graph G{Q~) (Corollary [38]). So the contracted 
subgraphs are precisely the 2-edge-connected components of G[Q. ■ 
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In summary we find all the V{B)-sets as follows. 

Compute all values ^„„, (^'^'" , uv an edge of G. Then repeat the following step for (^ taking on 

the distinct values in O in decreasing order: 

Construct G(C), contracting all V{B)-sets of blossoms of z- value > Q. Find 
the 2-edge-connected components of G{C,) and output them as the y(i?)-sets of 
blossoms of z- value C,. 

In addition, output the graph G{C,~) for the final value of C,. 

It is easy to modify the output to get most of the weighted blossom forest [{i)-{iv) of Seel?]) for 
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(i) The weighted blossom trees W are constructed from the containment relation for the 2-edge- 
connected components. 

(ii) z{V) is the final value of (. Consider a weighted blossom B ^ V. Let it be formed in the 
graph G{(). Then z{B) = C ~ C: where G(C') is the graph in which the parent of B (in W) is 
formed, or if B e M -V, (' = z{V). Clearly z{B) > 0. 

(in) The 7(i?)-sets are computed in the next section. 

(iv) The tree T is the final graph G{C~)- In proof G{C~) is a forest (Corollary [38]). Its vertices 
are the contractions of the blossoms of Ai, since every vertex in a critical graph is in a blossom 
(Lemma I33p . G{(^~) is a tree since a critical graph is connected. Finally every edge AB of G(C ) 
belongs to I (A) © I{B) since it is type 1 (and any blossom C has I{C) = Fc D S{C)). 

We conclude the section by estimating the time for this procedure. The values in are sorted 
into decreasing order in time 0(m log n). We use a set-merging algorithm to keep track of two 
partitions of V{G), namely the connected components of the G{C) graphs, and the 2-edge connected 
components. The total time for set merging is 0{m + v?). 

The total time for the rest of the procedure is 0{m+in?). In proof, there are < m iterations (|0| < 
m). An iteration that does not change either partition (because its new edges are contracted) uses 
0(1) time for each new edge. There are < 2n other iterations (since each of these iterations decreases 
the number of connected components or 2-edge connected components). Each such iteration uses 
linear time, i.e., 0(1) time per vertex or edge of G{C,). We complete the proof by showing that < n 
edges belong to > 1 graph G[C,). 

A type or 2 edge is in one graph G{C,). A type 1 edge first appears in G{C,) for ^ = min{^„, C,v}- 
Either it gets contracted in this iteration (and so is in just one G{C,) graph) or it is a 1(^ edge in any 
future G(C) graph that contains it. Any G{C,) has < n 1^ edges (Lemma I37p . 

8.2 Finding the /(B)-sets 

A (^-blossom is a blossom of W with z-value (". The iteration of the reduction for (" finds I{B) for 
each ("-blossom B. For any y G V{B) define the set 

I{y) = I{B)n6{y). 

So any blossom B has I{B) = IJ{-^(y) • V ^ ^(-^)} (this is part of the definition of blossoms |14| 
Definition 4.14]; it also easily follows from (iv) of the definition of the weighted blossom tree, since 
every maximum F'l^v respects every blossom). The iteration for ( computes the sets I{y) for all 
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vertices y with C^y = C,. Clearly we can combine these sets with sets I{B') known from previous 
iterations (for values C' > C) to find the /-set of each ^-blossom. 

We compute the /(y)-sets in two steps. For y and B as above (i.e., B a (^-blossom, y € V{B), 
(y = () define the set 

/o(y) = {xy : xy G 6{V{B)), (y > C" or (y = Cxy + Cx}- 

Lemma 40. /o(y) ^ ^{v)- 

Proof. Consider an edge xy G loiv)- Since xy G 6{V{B)) we have xy G I{y) iff xy G Fy {Fy respects 
B). The latter certainly holds if xy is type 2. So it suffices to show xy is not type 0, and xy G Fy 
if xy is type 1. 

Suppose xy is type 0. Lemma [35lfi) shows (^^ < <^^^, i.e., the first alternative in the set-former 
for /o(y) does not hold. If the second alternative holds we have (^ = (^y = Qxy, so Lemma [35lfi) 
implies x G V{B)^ contradicting xy G 5{y{B)). 

Suppose xy is type 1. The first alternative in the set-former for Io{y) implies Qy ^ C,^^ and the 
second alternative has C,x 7^ Cxy Both relations imply C,y = C,xy and xy G Fy (Lemma I35(iii)). ■ 

To find the remaining /-edges for G(C), define the set 

IE = {uv : uv joins distinct (^-blossoms, Q = Qu = Qv = Cmi}- 
Lemma 41. Any vertex y with Qy = Q has I{y) — loiv) ^ IE. 

Proof. Let B be the (^-blossom containing y. Suppose xy G /(?/)• Thus xy G Fy, making xy type 
1 or 2. If type 2, Lemma ISSTzz) shows (^^^ < (^.y. Furthermore with Qy = Q and x ^ V{B) it shows 
the inequality is strict. Thus xy G Io{y). If type 1, Lemma \35\ iii) shows (^y = Qxy If Cx 7^ Cy then 
xy G Io{y). If Cx = Cy then xy G IE, since xy G (5(y(/?)) shows x and y are in different (^-blossoms. 



IE is a subgraph of G(C~) (ti^ G IE has ^"^ = d^uv = C by (p3|) ). Now observe that every 
nonisolated vertex of G(C~) is a contracted blossom. (In proof, any edge xy of G{C~) belongs 
to G{C), so uim{(^xy,C^^} ^ C- Furthermore xy is type 1, so Lemma [35lfizi) shows min{(^a;,(^y} = 
in.m{(^xy, C^^}- Thus minj^^^, ^y} > C.) We conclude that every edge xy of G(C~) has xy G I{x)(Bl{y) 
(since every edge of G(C~) is type 1). 

So to complete the computation of the I{B)-sets we need only decide which alternative [xy G 
/(x) or xy G I{y)) holds for each edge xy of G{C~). (Note that an edge xy of G(C~) needn't be 
in IE - the blossoms containing x and y may not be (^-blossoms. But this causes no harm.) We 
accomplish this classification using the acyclicity of G(C~), as follows. 

Let T be a nontrivial tree of G(C~). Let i? be a leaf of T, incident to edge xy of T with y ^ B. 
Since Fy respects B, fy{B) + |/(/?)| is even (recall the definition of respect). Thus 
(*) xy G I{B) iff fy{B) + \I{B) - xy\ is odd. 

All edges of I{B) — xy are known. This follows from Lemma |4T] if ^ = (^y (recall I{B) — IE is 
known). If (^ < ^^ it holds since all of I{B) is actually known. So (*) can be used to add xy to 
exactly one of the sets /(x) or I{y). Thus the following algorithm correctly classifies each edge of 

Gin. 
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In the (current) forest G{C~), let i? be a leaf of a nontrivial tree T. Halt if no such T exists. 
Let xy, y £ B, he the unique edge of T incident to B. Assign xy to exactly one of I{x) or 
I(y), using (*). Then delete B from T and repeat. 

The total time to compute /(i?)-sets is 0{7n + n^). Specifically, a set Io{y) is computed in the 
iteration where Q = Qy, using 0(1) time on each edge incident to y. The algorithm for processing 
IE is executed in each iteration where some y has ^j^ = <^. Each such execution uses 0{n) time, 
giving time 0{n?') in total. 

8.3 Finding a maximum factor 

Assume we are given the dual functions y, and z in the form of its weighted blossom forest. We 
show how to find a maximum lower or upper perturbation F\:u for any given v £ V , va total time 
0{(f)^)^ with high probability. Let F be the desired maximum perturbation {F = F^ or F = F'"). 

The procedure is in three steps. It halts with the set Fq equal to F. Initially Fq is empty and 
we add edges to Fq as they become known. 

Edges of T 

The first step determines F OT and adds these edges to Fq. Consider any edge e of T- Let X and 
X' be the sets of the partition of V{G) induced by the connected components of T — e. Since F is 
an /;|:i;- factor, fv{X) counts every edge of F n "f{X) twice. So 

e G F iff f^{X) is odd. 
Clearly we can implement this test to find all edges of F Pi T in time Oiv?'). 

Strictly underrated edges 

The second step calculates yz{e) for each edge e G E{G). If e is strictly underrated it is added to 

To do the calculation efficiently assume each blossom is labelled with its z- value. Let e = uv. 
Let B be the nearest common ancestor of u and v in the blossom tree containing u and v. Then 

(ziB) uviI{Bu)yM{B,) 

yz{uv) = y{u) + y{v) + I z{Bu) uv G I{Bu) - I{B^) 

[ziBu) + z{Bv) - z{B) uv G I{Bu) n I{B^). 

This step uses total time 0{m). (Nearest common ancestors are found in 0(1) time.) 

It is convenient to ignore these strictly underrated edges in the rest of the discussion. So assume 

the degree-constraint function flp has been decreased to account for the strictly underrated edges, 

i.e., from now on f^p denotes the residual degree constraint. 

Edges of blossoms 

We turn to the third step of the procedure. As in shortest paths and 6-matching, we process the 
weighted blossoms in a top-down fashion. Consider a weighted blossom B. Assume the set Fr]6{B) 
is known. The first step ensures this for a root B £ Ai. We will find the edges of F that belong to 
G{B) and add them to Fq. Note that these edges complete the sets F n 6{A), A a child of B. So if 
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A is a blossom we can process it the same way. Thus we can process every weighted blossom this 
way. 

Form a graph H as G{B) with the strictly underrated edges deleted. Define e G E{G) U {0} by 
the relation 

Fr\5{B) =I{B)®e. 

li V ^ B this defines e as an edge, and if t; G i? it defines e as (recall F respects B). 

Note that if e is an edge it is tight. In proof, take any y G V{B). Then {Fr\5{B))®{Fyr]5{B)) = 
{I{B) ® e)) ® I{B) = e. So e G -F © -Fy implies e is underrated but not strictly underrated, whence 
e is tight. 

The desired subgraph of H is specified in Lemma [M] and we find it as follows. Let /' be the 
degree constraint function for H. Let x be a vertex of H. If x is a vertex of G then it has degree 
constraint 

f'{x) = fXv{x)-\5{x,Fr^5{B))\. 

If X is a contracted blossom A then the edges of F that are incident to x are governed by the 
fact that F respects A. Specifically F n 5{A) is I{A) if u G ^, and I{A) © g for some edge g G 5{A) 
\iv ^ A. 

|14| Section 4.4] models all the above constraints on the desired subgraph so that it corresponds 
to a maximum weight /-factor on H. It does this by redefining the edge weights (in fact it uses 
weights that are much larger than the given ones - see |141 Fig. 13]). This is inappropriate for the 
current context, since we wish to find the desired subgraph using a routine for unweighted /-factors. 

We model the constraint for A using the blossom substitute of FigJSj As indicated, the new 
vertices a,c and each a^ all have degree constraint 1, and 

, jo V G yl or e an edge in 5{A) 

1 1 otherwise. 

The following claim shows this substitute faithfully models the constraints on A. Let 5 be a set of 
edges in the substitute. Let the images of these edges in H be the set SB (e.g., ikdk in S corresponds 
to ikA in SB). Let 

SB = Fn S{B) n 6{A) = {I{B) © e) n 6{A). 

Claim S satisfies the degree constraints of the blossom substitute for A iff SB U SB respects A. 

Proof: Consider the two possible values of f'{b). 

Case f'(b) = 0: S satisfies the degree constraints iff it consists of edge ac and every edge a/^ii^ 
(but no edge bjk) i.e., SB = I (A) ^{B). This is equivalent to 

SBUSB = {I{A)®e)n6{A). 

li V £ A then e = 0, and the displayed equation becomes SB U SB = I{A) which is equivalent 
to SB U SB respecting A. Similarly \i v ^ A and e G 6{A) the displayed equation becomes 
SB U SB = I{A) © e which is equivalent to SB U SB respecting A. 

Case f'{b) = 1, v ^ A, and e ^ S{A): The assumption e ^ 6{A) implies SB = I{A) n d{B). 

The degree constraint /'(c) = 1 means S contains either edge ac or be. The former makes S 
contain each edge akik and one edge bjh, so SB U SB = I (A) + Ajh. The latter makes 5 contain 
one edge aoh and the edges a^-ik, k ^ h, so SB U SB = I{A) — Ai^. The assumption v ^ A in this 
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(a) i, 




(b) 




(a) Contracted blossom A in G{B). I and / (b) Substitute for A: vertices a, c, and each a^ all have /'-value 

denote 7(A) n 7(B) and 5(yl,i/) -J(yl) re- 1; & has /'(b) G {0, 1}. 

spectively. 

Figure 5: Blossom substitute. 



case shows that the two possibihties combined are equivalent to SB U SB respecting A. 







In summary the algorithm for the third step forms the above variant of H (note the tight edges 
have been identified in the second step). It finds an /'-factor on H and adds these edges to Fq. 

Lemma [Ml shows that the desired factor on H exists. (Note the hypothesis is satisfied, i.e., e is 
tight.) Now it is easy to see that the totality of edges added to Fq (by all three steps) achieves the 
optimality condition given in Section [3 for the desired maximum weight /;|:i)- factor. 

The time for the third step is 0((/)'^). To prove this it suffices to show that the total number 
of vertices in all H graphs is 0{(f)). There are 0{n) weighted blossoms, and hence 0{n) = 0{(f>) 
vertices of type a,b or c in blossom substitutes. An edge e in an /(i?)-set occurs in only one H 
graph - the graph corresponding to G(j>{B)), for B the maximal set with e G I{B) and p{B) the 
parent of B in W. e introduces 1 extra vertex a^ in H. A vertex u is on < /^(u) edges of sets I{B) 
(since these edges are in Fy) so the total number of edges e G [J{I{B) : B G W} is < (p. 

9 Shortest-path tree algorithms 

We construct the gsp-tree from the weighted blossom forest W found in Section [HI This section 
shows how to find a gsp-tree for the graph G{B), B a weighted blossom. It is a simple matter to 
join these gsp-trees together to get the entire gsp-tree for the given graph. (Note the root of the 
gsp-tree for any blossom B is a cycle node, not a tree node.) Let us restate the properties of G{B). 
Every edge is tight. There is a sink vertex t. Some vertices v have a distinguished edge e(f ) G 5{v), 
known to be on the shortest f t-path. (These vertices are children of B that are contracted blossoms, 
possibly singleton blossoms.) The remaining vertices of G{B) have no such edge (these vertices are 
children of B that are vertices in the given graph). We wish to find a gsp-tree for this graph G{B). 
We solve a slightly more general problem. Consider a graph with a sink vertex t and with every 
vertex v having a value eo{v) G 6{v) U {0}, eo(t) = 0. A -yu'-path P is permissible if every vertex 
X £ P — v' has eo{x) G E{P) U {0}. In particular P starts with eo(^) if it is nonnuU. Observe 
that any tit-path specified by a gsp-tree is permissible for e. We will present an algorithm that, for 
a given G, t,eo, finds a gsp-tree T whose function e agrees with eg on vertices where cq 7^ 0. In 
other words we will find a gsp-tree that specifies a permissible ui-path for every vertex v of G. For 
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Figure 6: Illustration for permissible paths: Arrows indicate e{v) edges. Solid edges belong to a 
search tree. The 3 leftmost dashed edges cause contractions; the 3 rightmost cause no action. 

convenience we assume the desired gsp-tree exists, as it does in our application. 

To illustrate the discussion Figl6] shows a given graph (all cq values are edges). Vertices i,j, k, m 
have no permissible path to t. FiglT] gives a gsp-tree for the subgraph of Figj6] induced by t and 
a, ... , h. 

We begin with some definitions and facts that lay the foundations of the algorithm. The discus- 
sion refers to both the given function eo and the final function e. In situations where either function 
can be used we try to use the more informative cq- Also note that the notion of permissibility is 
slightly different for the two functions. We will specify the function eg, e when it may not be clear. 

A cycle C is permissible (for e) if for some vertex b € C, every vertex x £ C — b has e{x) G E{C); 
furthermore, b = t if t £ C and b ^ t implies e{b) G S{C). Observe that for a permissible cycle C, 
any x £ C — b has a permissible path to b using edges of C. 

We will find the desired gsp-tree T by growing a "search tree", repeatedly contracting a permis- 
sible cycle C, and making C a cycle node of T. If C denotes the contracted vertex, the contracted 
graph has e(C) = e(6). Also if 6 = t then the contracted graph has t = C. 

As an example in FiglH] we will contract a, b, c; this gives the leftmost cycle node of FiglTj the 
contracted vertex gets e- value at. 

Let G be a graph formed from G by zero or more such contractions. A search irec£l T is a tree 
in G rooted at t, with each of its nodes v having one of two types, "u" for "up" or "d" for "down", 
depending on where we find the edge e{v). In precise terms let / be the edge from the parent of v 
to V. If / = e{v), or if f = t, then f is a u-vertex. In the remaining case (i.e., v ^ t and / ^ e{v)) v 



The reader will recognize the resemblance to search trees used in cardinality matching algorithms. 
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Figure 7: Gsp-tree for vertices a, . 
children of A^. 



, h and t of Figl6l e(A^) edges are the dashed edges joining the 



is a d-vertex. We require that a d-vertex v has at most one child, namely c where e{v) = vc. 

As an example in FigiG] a search tree might initially contain u- vertices t, a and d- vertices b, c. 
Contracting a, b, c makes the contracted vertex a u- vertex. Edge cf can be added to the search tree 
after this contraction, but not before. Adding cf makes / a d-vertex. 

A search tree has the following property: 

Lemma 42. Consider a nontree edge X1X2 where for i = 1,2, Xi is a u-vertex or Xi is a d-vertex 
having e(xj) = 3;i2;2. Let b be the nearest common ancestor of xi,X2 and let C be the fundamental 
cycle of xiX2- Then C,b is a permissible cycle. 

FiglH] illustrates the lemma. In the given graph (i.e., before any cycles have been contracted) 
edges be and et give permissible cycles (for an appropriate search tree), be (et) illustrates the lemma 
when {3;i,3;2} contains 2 (1) d- vertices, respectively. After these 2 permissible cycles have been 
contracted, gd illustrates the lemma when {xi,X2} contains 2 u- vertices. The nontree edge ik (jl) 
illustrates a nonpermissible cycle, when {xi,X2} contains 2 (1) d- vertices, respectively. 

Proof. First note 6 is a u-vertex. (If 6 is a d-vertex, it has just one child c, with e(6) = be. 2;i2;2 
must be a back edge with some Xj = b. The hypothesis makes e(xj) = e(6) = X1X2, contradiction.) 
Now it is easy to check the properties for permissibility (e.g., a d-vertex v £ C has e{v) € C since 
V is either an Xj or the edge from v to its unique child is in C). ■ 

The algorithm will contract such cycles C, and only these. In particular this implies that in G, 
any d-vertex is actually a vertex of G. 

Also note the algorithm will choose values e{v) for vertices with eo(f) = by making v a 
u-vertex. In particular this implies that e(i') = eo(i') for every d-vertex v. 
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Algorithm We will build the desired gsp-tree by growing a search tree T. Initialize T to consist 
of just a root t (a li-vertex) and let e be the given function cq. Then repeat Algorithm [3] until every 
u- vertex x in the current graph has every incident edge either in T or scanned from x, and every 
d-vertex has been scanned. 

Algorithm 3 A procedure that scans vertices and builds the gsp-tree. 

1: if some u- vertex x has an unscanned edge xy ^ T then > scan xy from x 

2: if y ^ T then 

3: add xy to T 

4: else if y is a u- vertex or xy = eo{y) then \> y is a d-vertex leaf 

5: contract the fundamental cycle of xy 

6: end if 

7: else if there is an unscanned d-vertex x then > scan x 

8: let 60(3;) = xy 

9: if y ^ r then 
10: add xy to T 

11: else if xy = eo(y) then \> y is a d-vertex leaf 

12: contract the fundamental cycle of xy 

13: end if 
14: end if 

Addendum Several aspects of this algorithm are stated at a high level and deserve further eluci- 
dation: 

An edge xy that is scanned from a u-vertex x is considered scanned from any u-vertex x that 
contains x by contractions. 

As mentioned, when the algorithm adds an edge xy to T, if eo(y) = it sets e{y) = xy. 
Furthermore in all cases it makes y type u or d as appropriate. 

When the algorithm contracts a fundamental cycle C, it creates a node in T whose children 
correspond to the vertices of C If the algorithm halts with T containing nodes other than t, it 
creates a root node of T whose children correspond to the vertices of T. 

Examples Note that over the entire algorithm an edge may get scanned twice, once from each 
end, e.g., in Figl6] (7(i may get scanned from g when d is still a d-vertex. In a similar vein, if we 
change eo{i) to a new edge id, the algorithm might scan i (and edge id) before d becomes a u-vertex 
and id gets scanned from it. 

Analysis of the algorithm The algorithm takes no action for the following types of edges xy: 

(i) X a u-vertex, y a d-vertex, xy 7^ eo(y), 

(ii) X and y d- vertices, xy is not both eo{x) and eo{y). 

(Hi) X a d-vertex, y ^ T, xy ^ eo{x). 
(ii) and (Hi) include edges that are never even scanned (e.g., ik,im in FigiG]). When the algorithm 
halts every nontree edge of G with an end in T is of type (i), (ii) or {Hi). 

Now we show the algorithm is correct. First note Lemma |32] shows every contracted cycle is 
permissible. 
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Assume every vertex v has a permissible path to t in the given graph G with function cq- Let 
G be the final graph of the algorithm. Let T be the final tree. 

Lemma 43. G consists entirely of u-vertices. 

Proof. Let v be a non-u-vertex, with P a permissible tit-path in G. The reader should bear in mind 
the possibility that the image of P in G needn't be permissible (for e). For instance suppose G 
contains a contracted node like X = {a, ... , g, t} in FigjGj If the graph has two edges b'b, c'c € 5(X), 
P might contain a subpath b',b,c,c', making it nonpermissible in G. In fact similar edges incident 
to fg might make the image of P nonsimple in G. 

First observe that we can assume w is a d- vertex and V{P) C y(T)l^ In proof, suppose P 
contains a nontree vertex (if not we're clearly done). It is eventually followed by an edge rs, where 
r ^ V{T) and V{Q) ^ ViT) for Q the st-subpath of P. {i)-{iii) show s is not a u-vertex. So s is a 
d- vertex. [Hi) shows rs ^ eo(s), so Q is permissible. 

Imagine traversing the edges of P, starting from v. Some edges of P will be in contracted 
vertices of G, others will be edges of T, and all others will be of type {i) or [ii), not [Hi) (since 
V{P) C V{T)). We assert that whenever we reach a vertex r of G, either 

(a) r is a d- vertex, and eo(r) has not been traversed, or 

(b) r is a u-vertex, reached by traversing edge e(r). 

We prove the assertion by induction. Note that the assertion completes the proof of the lemma: 
(b) shows P always enters a u-vertex from its parent, so it never reaches the u-vertex t, contradiction. 

For the base case of the induction, r = v obviously satisfies (a). For the inductive step assume 
the assertion holds for r and let s be the next vertex of G that is reached. So rs is the next edge 
of G that is traversed, (r or s may be contracted vertices of G.) 

Case r satisfies (a): Permissibility in G implies the next edge of P in G is rs = eo(r). 

Suppose eo(r) G T. This implies s is the child of r. The definition of search tree implies s 
satisfies (a) or (b), depending on whether or not rs = eo(s). 

Suppose eo{r) ^ T. eo{r) is not type (i) above (even if we take r = y). So it is type [ii). This 
implies (a) holds for s. 

Case r satisfies (h): The argument is similar. 

Suppose rs G T. (b) shows r is reached from its parent in T. So s is a child of r. As in the 
previous case, s satisfies (a) or (b). 

Suppose rs ^T. So rs is type (i) above and (a) holds for s. ■ 

The lemma implies T with its contractions gives the desired gsp-tree T. (E.g., if T consists of 
just one vertex t, the root of T is a cycle node, specifically the last cycle to be contracted; otherwise 
the root is a tree node and its tree is T.) 

It is easy to implement the above procedure in time 0{m\ogn) or better using an algorithm for 
set merging |11) . 



^^Here we commit a slight abuse of notation: The inclusion is meant to allow the possibility that a vertex of P is 
contained in a contracted vertex of T. 
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10 Combinatoric algorithms for shortest paths 

Let (G, t, w) denote a connected undirected graph with a distinguished vertex t and a conservative 
edge- weight function w : E ^ M. Let E^ be the set of edges with negative weights. In this section 
we show how to use combinatoric algorithms for finding maximum perfect matchings to compute 
the gsp-tree. We will define a graph Gt that models paths in G by almost perfect matchings, i.e., 
matchings that miss exactly one vertex. Moreover, we define a v-m,atching to be an almost perfect 
matching in G that avoids v. We believe that the construction of the split graph G is essentially 
due to Edmonds [9]. We define the split graph G = {V,E) with weight function id in the following 
way 

y = {7;i,t;2 -v G V}U {61,62 : e G E~}, 

E = {f 1^2 : V £ V} U {uiV2,U2Vi,uiVi,U2V2 : uv £ E\ E~} 

U {ui6l,U26l,6l62,Vl62,V2e2 : 6 = UV £ E~ , U < v} , 



w{UiVj 



-w{uv) if uv £ E\ E~ , 
—w{e) if Ui = ei and Vj ^ 62 and 6 £ E' 
otherwise. 



An important property is that we can assume fi = |1/| < 4n. This follows since we can assume 
\E~\ < n, as otherwise the set of negative edges contains a cycle. The following observation is 
essentially given in |T] in Chapter 12.7, where the reduction is explained on a clear example. 

Lemma 44. Let u,v £ V, let M be the maximum perfect matching, and let M{u2Vi) he the max- 
im,um weight perfect matching in G — U2 — vi. If G does not contain negative weight cycles then 
'w{M) = and the shortest path weight from u to v in G is equal to —w{M[u2Vi)). 

Note also that it is easy to detect a negative cycle in G - it corresponds to a perfect matching in 
G with positive weight. On the other hand, as described in [1], in order to find a shortest path from 
ti to i; we need to find the maximum weight perfect matching M{u2Vi). However, here we want 
to compute the whole gsp-tree and hence require the distances from all vertices in G to t. We will 
show that in order to find all these distances it is essentially enough to find one maximum perfect 
matching. Let us define Gt to be graph G with both vertices ti and t2 unified to one vertex t. We 
observe that the resulting graph is critical, i.e., for each vertex v there exists a i;-matching. 

Lemma 45. Graph Gt is critical. Moreover, let M{vi) he the maximum vi-matching in Gt then 

w{M{t2Vi)) = w{M{vi)). 

Proof. We need to show the existence of u-matchings for all vertices f in Gf . Consider the following 
cases 

• V = t then M{t2ti) in G corresponds matching that avoids t, 

• V = V2 then M[tiV2) in G corresponds to matching that avoids V2, 

• V = vi then in this case we take M(t2^i) in G. 

Hence, Gt is critical. The second part of the lemma follows by the above correspondence of match- 
ings. I 
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We need to relate the above observation to the definitions from previous sections. You might 
observe that M{vi) corresponds to Py from Section [6] and in language of /-factors to Fy as discussed 
in Section [8j In order to construct the gsp-tree we use the shrinking procedure for general /-factors 
from Section [9l This requires us to know w{F^) as well. However, we can observe that F'" can be 
obtained from Fy and vice versa by adding or removing the zero weight loop vv. 

10.1 Matching duals for factor critical graphs 

Consider an arbitrary graph G = (V, E) that contains a perfect matching. Let zz; : i? — )• M be 
the edge weight function. The dual variables in Edmonds' formulation [8] are assigned to vertices 
y : y — 7- M and to odd-size subsets of vertices 2; : 2^ — )• M. The function z can be negative possibly 
only on V . We define the value of uw (with respect to the dual y, z) as 

yz{uv) = y{u) + y{v) + z{B : e C B}. 
We require the duals to dominate all edges uv G E, i.e., we require 

yz{uv) > w{uv). 

We say that an edge is tight when the above inequality is satisfied with equality. The dual objective 
is defined as 

(y, z)V = y{V) + Y,{WB\/2\z{B) : B C F} . 

By the duality we know that for any perfect matching M and duals y, z we have w{M) < (y, z)V. We 
say that a matching respects a set B if it contains [|i?|/2j edges in "){B). As shown by Edmonds J8] 
a perfect matching is maximum if and only if all its edges are tight and it respects all sets with 
positive z for a pair of dominating duals. 

A blossom is a subgraph B oi G defined as follows. Every vertex is a blossom and has no edges. 
Otherwise, the vertices V{B) are partitioned into an odd number k of sets V{Bi), for 1 < i < A;, 
where each Bi is a blossom. Each blossom B contains edges in E(Bi) and k edges that form a cycle 
on Bi, i.e., consecutive edges end in V{Bi) and V{Bi^i), where Bi = fij+i. Blossoms Bi are called 
subblossoms of B. The set of blossoms can be represented as forrest called blossom, forrest. The 
parent-child relation in this forrest is defined by the blossom-subblossom relation. 

Edmonds' algorithm for finding maximum perfect matchings constructs a structured matching, 
i.e., a matching M and a dual solution. The dual solution is composed out of a blossom forrest F 
and functions y and z. The structured matching satisfies the following conditions 

(i) M respects blossoms in F, 

(ii) z is nonzero on blossoms in F, 

(Hi) all edges in M are tight, 

(iv) all edges in blossoms in F are tight. 

If one finds a structured matching that is perfect then it is a maximum perfect matchings. In 
the classical view Edmonds' algorithms algorithm operates on a graph with even number of vertices. 
However, as shown in |15| it can be seen to work on a critical graph. The optimal matching structure 
of a critical graph G consists a blossom tree B and dual functions y, z such that every vertex is 
a leaf in B and properties (ii) and (iv) above are satisfied. Observe that there is no matching in 
this definition. As argued in |15| Edmonds' algorithm computes an optimum matching structure 
when it is executed on a critical graph. Moreover, the optimal matching structure allow us to relate 
weights of maximum f-matchings to the dual y, z in the following way. 
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Lemma 46 (|15|). Let y,z and B the optimal matching structure for the critical graph G. Then 
the weight of maximum v-matching is to equal to {y,z)V — y{v). 

10.2 The algorithm 

Let us now join all the ingredients to compute the gsp-tree in the undirected graph (G, w, t) with 
the conservative weight function w. 

Algorithm 4 A coinbinatoric algorithm for computing gsp-tree for the undirected graph {G,w,t). 

1: Construct G from G 

2: Construct Gt from G by identifying ti and ^2 

3: Compute optimal matching structure y,z and B for Gt > Gt is critical by Lemma l45l 

4: for V E F do 

5: Let w{My) = (y, z)V — y{v) t> by Lemma B6l 

6: Let w{F^) = —w{M^) > by Lemma HSJ and Lemma |4H 

7: Let w{F^) = w{Fy) l> by adding zero length loop vv 

8: end for 

9: Using w{Fy) and w{F^') find the blossom forrest > using shrinking procedure from Section [8] 
10: Construct gsp-tree from the blossom forrest > using procedure from Section [9] 

We note that optimal matching structure can be found using fast implementations of Edmonds' 
algorithm in 0{n{m + nlogn)) time |13| or in 0{y^na{m, n) log n mlog{nW)) time |15| . The 
other steps of the above algorithm take only less time, e.g., the shrinking procedure can be im- 
plemented to work in 0(m log n) time. This way we obtain the 0{n{m + nlogn)) time and the 
0{y^na{m, n) log n mlog {nW)) time combinatoric algorithms for computing gsp-tree. 

11 Determinant formulations for general graphs 

Let G be a simple graph with vertices numbered from 1 to n. Let cj) = ^^ /(*)• We define a skew- 
symmetric (p X (f) matrix B{G) representing G in the following way. A vertex i ^ V \s associated 
with f{i) consecutive rows and columns of B, both indexed by the pairs i,r for < r < f{i). Call 
a tuple {i,r,j,c) representing the entry B{G)i^r,j,c permissible if either 

i < j, or i = j and r < lf{i)/2\ < \f{i)/2] < c < f{i). 

The set of permissible entries is denoted by P. Note that the permissible entries of B[G) are all 
above the diagonal. The permissible entries corresponding to a fixed edge e form a rectangular 
submatrix for a nonloop e and an \_f{i)/2\ x [/(i)/2j submatrix for a loop ii. 
Using indeterminates x'r , Vc we define an entry of B{G) as 



B{G) 



i,r,],c 



XrVc 


ij € E and (i, r, j, c) € 


-x'cy'r 


ij G E and (j, c,i,r) £ 





otherwise. 



(25) 



Clearly B{G) is skew-symmetric. Note that a loop ii £ E is represented by an [/(i)/2j x [f{i)/2\ 
submatrix that is empty if f{i) = 1. This is fine since ii is not in any /-factor. 
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The analysis for multigraphs is almost identical to simple graphs. So we will concentrate on 
simple graphs, but also point out how it extends to multigraphs. Towards this end we extend the 
definition of B{G) as in Section [4.11 Let //(e) denote the multiplicity of any edge e. The copies of 
e = ij get indeterminates x'r' , 



yT'' {l<k< //(e)) and we define B{G) by 




' Efcil xV'^Vc'^ ij G E and (i, r, j, c) G 


< 


-ES^PyP iiGii;and(j,c,/,r)G 




otherwise. 



B{G\r,,c = { - T.tn ^'c'y'r''" ij G E and (j, c, z, r) G P, (26) 



11.1 Review of the PfafRan 

Our goal is to prove an analog of Theorem [5l It is easiest to accomplish this using the Pfaffian, 
so we begin by reviewing this concept. Let ^ be a skew-symmetric matrix of order 2h x 2h. Its 
Pfaffian is defined by 



Pf (^) = X] '^Sn 



1 2 ... 2^-1 2h 
ii ji ■■■ ih jh 



'■iiji ■•• "-ihjh- 



The above sum is over all partitions of the integers [1..2/i] into pairs denoted as {ii,ji}, ... .{ih,jh}- 
Also sgn denotes the sign of the permutation. (It is easy to check that each term a in the summation 
is well-defined. For example interchanging ii and ji does not change the partition. It fiips both 
the sign of the permutation and the sign of aj^ji, so a is unchanged.) The central property of the 
Pfaffian is [H] 

det(A) = {pi{A)f. (27) 

Recall that the Tutte matrix T for a graph G is the skew-symmetric matrix obtained from the 
adjacency matrix of G by replacing the entry 1 for edge ij by an indeterminate tij if i < j and by 
—tji ii i > j. Clearly there is a 1-1 correspondence between perfect matchings of G and terms of 

pf(r). 

Recalling ()27p . note that every pair of perfect matchings of G gets combined in the product 
(pf (T))^. Let us review a proof that no cancellations occur when all these terms are added together. 

Consider a term of (pf (T))^ corresponding to perfect matchings Mi, M2. (These matchings may 
be distinct or identical.) The terms of (pf(r))^ that might cancel this term must use exactly the 
same variables tij . In other words they correspond to perfect matchings Ni , A'^2 where 

Ml tt) M2 = iVi l±) N2. (28) 

Here tt) denotes multiset sum. Mi tt) M2 consists of the edges of Mi n M2 taken twice, plus the set 
Ml M2. The latter is a collection C of even alternating cycles. The definition of matching shows 
every cycle of C is alternating wrt Ni,N2- Assume that in (pf(T))^, the first multiplicand pf(T') 
gives the matching A'^i and the second multiplicand gives A'^2 ■ The edges of each cycle of C can be 
partitioned in two ways between A'^i and A'^2- So there are precisely 2' ' matching pairs Ni,N2 of 
(pf(r))2 that satisfy (l28]l . 

We will show there are no cancellations because each matching pair Ni,N2 gives the same term 
of (pf(r))^. More precisely we show the following: 

Claim 1. The pairs Ni,N2 satisfying (|28|) collectively contribute the quantity (— 2)' 'lltjj to 
(pf(T))^, where the product is over a fixed set of entries in T corresponding to Mi tt) M2. 
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Note the claim also implies that if we do arithmetic over a finite field of characteristic > 2, again 
there are no cancellations. 

Proof of Claim 1. Take a pair of matchings Ni,N2 satisfying (j28p . Let CTj be the term in pf(T) for 
Ni, so (pf(T))^ contains 0"i(72. We will show the portion of 0"i<T2 corresponding to C is the same 
for every Ni,N2. Then we will conclude this property makes the entire term (Ti(T2 independent of 
choice of Ni,N2. 

First suppose C is a single edge ij (belonging to Mi M2). Wlog both o"! and (T2 have the T 
entry tij and their permutations both map some ordered pair of two consecutive integers (2a— 1, 2a) 
to the ordered pair {i,j). Clearly this does not depend on choice of iVi, N2- 

Now suppose C contains vertices ii, ... , ii for some even ^ > 4, and A'^i contains edges iii2, ••• , h-iH 
while A'^2 contains ^2^35 ••• ) HH- <7i contains the product ij^jj ■■■iie-i,ie ^-i^d o"2 contains tjjjs ■■■tie,ii- 
So o"i(T2 contains the product {ti^i^ ■■■tii_i,ii){ti2i3 ■■■tie,ii)- Certainly this is independent of choice 
ofiVi,iV2. 

The permutation for ui maps some pair of two consecutive integers (2a — 1,2a) to (ii,i2)) and 
similarly for the rest, e.g., (26 — 1,26) goes to (ii-i,ii). Wlog A'^2 maps (2a — 1,2a) to (i2i^3)i 
and similarly for the rest, e.g., (26 — 1,26) goes to (if,ii). The pairs of cJi are transformed to the 
pairs of (T2 by applying the product of transpositions (1112) {iiis) ■■■ (hie)- This is an odd number of 
transpositions, i.e., its sign is —1. So the part of the permutations for o"! and o"2 hi C combine to 
give the sign —1 in 0"icr2. 

Applying this analysis of sign to each of the \C\ components of > 4 vertices shows the entire 
permutation of (Ji(T2 contributes sign {—!)''. Since there are 2' ' pairs Ni,N2, we get the contri- 
bution of the claim. (} 

We shall also use a special case of the above analysis, specifically when A'"i and N2 are identical 
except for choosing alternate edges of one cycle of length 4. ui and (T2 have opposite sign, and so 
they differ only in the subexpressions 

for some consistent choice of sign. 

11.2 Analysis of B(G) 

We return to the matrix B{G) for /-factors. An entry of B{G) corresponds to an edge of G, so a 
term of pf (i?(G)) corresponds to a multiset of edges of G. 

Claim 2. Each uncancelled term a in pf(i3(G)) corresponds to an f -factor of G. 

Proof of Claim 2. The above discussion shows that ignoring sign, cr is a product of (j)/2 quantities 
Xr^yl^ corresponding to edges that form a matching on vertices designated by two indices i,r (1 < 
i < n, < r < f{i))- So each vertex i of G is on exactly f{i) of the (p/2 edges selected by a 
(counting loops ii twice). To show these edges form an /-factor of G we must show that a uses 
each edge of G at most once. We accomplish this by showing that in the summation of the Pfaffian, 
terms using an edge more than once cancel in pairs. 

Suppose a uses edges Xryl^ and x^^,y^^, (r / r' and c^^ d). Wlog assume (i, r; j, c), (i, r'; j, c') G P. 
a may have many such duplicated pairs. Choose the duplicated pair that lexically minimizes (i, r, r'). 
Pair a with the term a' having the same partition except that it replaces {i, r; j, c} and {i, r'; j, c'} by 
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{i,r;j,c'} and {i,r';j,c}. Note that permissibility of {i,r;j,c) and {i,r';j,c') implies permissibility 
of {i,r;j,c') and {i,r';j,c), even if i = j. Thus a' is also associated with {i,r,r'), so the pairing is 
well-defined. 

Viewed as matchings, a and a' differ only by choosing alternate edges of the length 4 cycle 
{i,r;j,c;i,r';j,c'). So a and a' are identical except for the expressions of (j29p . which in the new 
setting become 

^"i,r;j,c"i,r';j,c' and -{- "j,c]i,r' ^j,c';i,rj 

where b designates matrix B[G). Substituting the definition of B{G) shows these expressions are 

±{+a^yi^){+:^,y>^) and ^ (-x>^)(-x^2/^). 

The two minus signs in the second expression follow from permissibility of (i, r; j, c') and (i, r'; j, c). 
The 4 indeterminates in the 2 above expressions are collectively identical. So we get a = — o"', i.e., 
these two terms cancel each other as desired. 

This argument also applies to multigraphs G. Here the issue is that only one copy of each 
distinct edge ij, k can used. A term a with a duplicated edge, like xl^' yl^' and x^^,' y^^,' , chooses 
the duplicated pair to lexically minimize {i,k,r,r'). The rest of the argument is unchanged. 

Every /-factor F oi G has an uncancelled term a in pf(i3(G)): For each z G y, order the set 
6{i,F) arbitrarily. This makes F correspond to a term a in pf(i?(G)). Each variable xl^ , y^c in cr 
specifies its edge, and so determines a unique partition pair {z, r; j, c}. So no other term of pf (i?(G)) 
has the same variables of u, and a is uncancelled, (Of course a given F gives rise to many different 
uncancelled terms.) 

Now it is easy to see that (f27|) gives an analog of the second assertion of Theorem [5) G has an 
/-factor if and only if det{B{G)) ^ 0. But we need to prove the stronger first assertion. 

We extend the above construction for F to an arbitrary pair of /-factors Fi,F2. As before for 
each i £ V and for each Fj (j = 1,2), number the edges of 6{i, Fj) (consecutively, starting at 1). But 
now for both sets start with the edges of 5{i,Fi D F2), using the same numbering for both. Using 
both these numberings gives two terms in pf (i?(G)), which combine to give a term a of (pf (i?(G)))^. 
Each variable x^r in o" specifies its edge and corresponds to a unique y^c . Here we are using the fact 
that an edge in both /-factors gives rise to a subexpression of the form (xl^yl^)'^ in a. 

Call a term of (pf(i?(G)))^ consistent if whenever an edge ij occurs twice, it uses the same 
subscripts, thus giving a product of the form {xryl'')'^. (We are disqualifying terms with an edge 
appearing as both xl^yc and also x''^,y^J, where r 7^ r' or c 7^ c'.) Clearly the above construction 
gives consistent terms. We will show there is no cancellation in (pf(i?(G)))^ involving consistent 
terms. (After that we comment on inconsistent terms.) 

The proof follows the Tutte matrix argument: Consider a pair of /-factors Fi,F2, and some 
consistent term a in (pf(i?(G)))^ that corresponds to them. View these /-factors as matchings 
FMi,FM2 on vertices i,r. Another pair of /-factors Hi^ H2 involves the same set of indeterminates 
exactly when FMi tt) FM2 = HMi tt) HM2. (Here we use the fact that by consistency, each 
indeterminate x^ specifies its edge ij as well as its corresponding indeterminate y^c . Hence the 
indeterminates in a term determine the matched edges (z,r;j, c).) We have now established the 
analog of (I28p . The rest of the argument for Tutte matrices applies unchanged. Furthermore 
exactly the same analysis applies when G is a multigraph. 

We comment that the situation is more involved for an inconsistent term a. Suppose Fi has 
Pfaffian term cJi containing x^r yi and F2 has G2 containing x^,J,y^i ^ where r ^ r' and c ^ d . 
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Furthermore suppose the two copies of ij belong to different cycles Ci, C2 of FMi tt) FM2- There 
can be many different possibilities for two other /-factors Hi,H2 having HMi tt) HM2. For instance 
suppose edge ij gives the only inconsistency in (Tia"2. Then Ci and C2 only give rise to 2 possible 
matchings HMi,HM2 instead of 4, since an /-factor cannot contain both copies of ij. But if there 
are inconsistent edges besides ij, other partitions for the edges of Ci and C2 may be possible (e.g., 
consider the case of 1 other inconsistent edge, with copies in Ci and C2). 

As in Section m define F to be the function that maps each term a of det{B{G)) to its corre- 
sponding subgraph denoted F^-. F„ is a 2/-factor (unlike Section H]). Let $2(G) be the set formed 
by taking sums of two (possibly equal) /-factors of G, i.e., 

$2(G) = {Fi W F2 : Fi, F2 /-factors of G}. 

We have proved the following. 

Theorem 47. Let G be a simple graph or a multigraph. The function T from terms in det(i?(G)) 
is a surjection onto ^2{G)- Consequently, G has an f -factor if and only ifdet{B{G)) 7^ 0. 

The theorem continues to hold when we do arithmetic in any finite field of characteristic > 2. 

12 Finding /-factors in general graphs 

This section gives algorithms to find /-factors in general multigraphs. It starts with simple graphs 
and then moves to multigraphs. One should keep in mind that it is unknown how to use Gaussian 
elimination in the non-bipartite case, e.g., [21] uses a different non-algebraic algorithm for this 
case. However Harvey |17| later developed a fully algebraic scheme, that we adopt here. A good 
explanation of this approach is given in [17]. We first use the Sherman- Morrison formula to get an 
OinKfP') time algorithm. Then we show it can be sped up to 0{(j)^) time using ideas from Harvey's 
recursive elimination scheme. 

12.1 Simple Graphs 

We define removable edge ij G E[G) to be an edge such that G — ij has an /-factor. We can observe 
that following property. 

Corollary 48. Let G be a simple graph having an f -factor. The edge ij £ E[G) is removable if 
and only if det {B{G — ij)) 7^ 0. 

Let ij G E{G) then we can observe that B{G) and B[G — ij) differ from one another by two 
rank-one updates. 

B{G - ij) = B{G) - x'\y'^)^ + y'^{x'^)'^ = B{G) + [-x'^ ,y%y'^ , x'^f. (30) 

where x^^ and j/*-' are length (f) vectors. When i ^ j these vectors are defined as 

xl^ \ik = j, j,- J yl!' iik = j, 



1 otherwise, '"^ \ otherwise. 
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On the other hand for ii G E{G) we have 

jj _ jxf if k = i and r < [f{'i)/2\ , ■■ _ J y** ii k = i and f{j)/2 < r, 

1 otherwise, '^ 1 otherwise. 

Hence, one can use Sherman-Morrison- Woodbury formula to compute B{G — ij)~^ from B{G)~^ 
in 0{(jP') time. Similarly, we can use this formula to test whether B{G—ij) is nonsingular by checking 

det(l2 + [y'^,x'^fB{Gr^ [-x^^y'^■]) / 0. (31) 

Here, the affected matrix size is bounded by 2f{i) + 2/(j) so we need 0{{f{i) + f{j))'^) time for 
this test. Using these observations we get Algorithm [5] for finding /-factors in simple graphs. 

Algorithm 5 An 0{m(fP') time algorithm for finding /-factor in the simple graph G. 

1: Let B[G) he (p X (j) skew-symmetric adjacency matrix of G 

2: Replace the variables in B{G) for random elements from Zp for prime p = @{(f)'^) to obtain B 

3: If B is singular return "no /-factor". 

4: (with probability > 1 — 4 matrix B is non-singular when B{G) is non-singular) > by Lemma [T] 

5: Compute B~^ 

6: (we remove only removable edges so B remains non-singular during execution of the algorithm) 

7: for all ij ^ E do 

8: if det(/2 -I- [y*-', x*-'] B~^[-x^^ ,y'^^]) ^Q then t> Edge ii is removable bv (f3T]) 

9: Set E := E - e 

10: Set B := B+ [-x'^ ,y'^] ■ [y'^,x'^f > This corresponds to B{G) - ij by ^ 

11: Recompute B~^ > Using Sherman-Morrison- Woodbury formula 

12: end if 
13: end for 
14: Return E > All removable edges have been removed, so what remains is an /-factor 

Observe that the above algorithm fits into the framework introduced by Harvey |17] for find- 
ing 1-factors in general graphs. The above algorithm corresponds to the algorithm described in 
Section 3.3 from |17| . with the difference that the rank-two updates we use affect submatrices and 
not single elements. This means that we can use the recursive algorithm that was introduced in 
Section 3.4 of his paper. The algorithm uses following elimination procedures with the starting call 
to DeleteEdgesWithin(y). 

DeleteEdgesWithin(S') - if \S\ > 1 split S into Si and ^2; call DeleteEdgesWithin(S'i), for i = 1,2; 
call DeleteEdgesCrossing(S'i, 52); update submatrix B~^[S,S]; 

DeleteEdgesCrossing(i?, S) - if \R\ = {r}, IS"! = {s} and rs is removable eliminate edge rs; otherwise 
split R into Ri,R2 and S into Si, S2', call DeleteEdgesCrossing(i?j, S'j), for i,j = 1,2; update 
submatrix B-^[R U S,RUS]. 



Observe that we need only to replace tests for removable edges with (I3ip . whereas the updates to 
submatrices remain essentially the same and use |17| Corollary 2.1]. Hence, the submatrix of size 
\S\ X \S\ is updated in 0{\f{S)\'^) time. On the other hand, the cost we pay to test whether the 
edge ij is removable is 0(|/(i) -|- /(j)P). 
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Harvey splits the set S (similarly R) always in equal halves. We, however, split S in such a way 
that f{Si) and f{S2) are as close as possible. Let us assume f{Si) > f{S2) then 

• either f{Si) > |/(5) and |Si| = 1, 

• or /(5i) < |/(5) and /(Sa) > i/(5). 

Let h{S) denote the running time of the procedure DeleteEdgesWithin(5'). Similarly define g{R, S) 
for DeleteEdgesCrossing(i?, 5). We have 

h{S) = ^ h{S,) + g{Si,S2) + 0{\f{S)n 

i 

.^5) = /0(l/(^) + /(5)P) if|/?| = |S| = l, 

' \Ei,j9{Ri,Sj) + 0{\f{R) + f{S)n otherwise. 

The solution for these equations gives an 0{(j)^) time bound. 

12.2 Multigraphs 

In the case of multigraphs we need to handle multiple copies of the same edge in a different way. 
Removing separate copies of an edge one by one would lead to a cubic time complexity. The cost 
charged by edge ij to its submatrix would be Q,({f{i) + f{j))'^ • fi{ij)). Instead, we use binary search 
on the number of removable copies of edges. Assume that we want to remove /_f copies of edge 
ij G E[G) from the graph. We denote the resulting graph by B{G — ij^)- In such case B{G) and 
B[G — ij^) differ from one another by 2fi updates of rank one, i.e., 

B{G - if) = B{G) + V -x*^>(2/*^''^)^ + y^^'>^{x'^''^f 

k=i ^'^^^ 

= B{G) + [-x*J'\ 7/^^'\ ... , -3;^J>, 2/*J>][y^J'\ x^J'S ... , y^^>, x^^>]^. 

where x*-'' and y^^' are length (j) vectors defined in similar way as in the previous section. Observe 
that the submatrix affected by these changes has size 2/(z) + 2/(j) so, computing B{G — ij^)~^ 
from B{G)~^ can be realized using [171 Corollary 2.1]. Moreover, we can test non-singularity of 
B{G- if) by 

det {I, + [-x'^^\ ... ,y*^'^]^i?(G)-My*^'\ ... ,x^^'1) / 0. (33) 

This test requires 0{[f[i) + f{j))'^) time, because the affected matrix size is 2/(i) -|- 2f{j). In the 
following algorithm we use this test together with a version of binary search. 
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Algorithm 6 An 0{m<j)^) time algorithm for finding /-factor in the multigraph G. 



2: 
3: 

4: 

5: 

6: 

7: 

8: 

9: 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 



Let B{G) he (p X (j) skew-symmetric adjacency matrix of G 

Replace the variables in B{G) for random elements from Zp for prime p = 0((/>^) to obtain B 

(with probability > 1 — 4 matrix B is non-singular when B{G) is non-singular) > by Lemma [T] 

Compute B~^ 

(we remove only removable edges so B remains non-singular during execution of the algorithm) 

for all ij € -E do 

Let /x be the highest power of 2 not higher then min(/i(ij), /(i), f{j))- 
Let k := 0. 
v^hile /x > 1 do 

if det {If, + [-x^J'^'+i, ... ,y*J>+^]^5-^[2/*J>+i, ... ,x*J>+/^]) / then 

> ^ copies of ij are removable by (|33p 
Set E:=E- if 
Set B := B+ [-x'^'''+'^ ... ,y*J>+/^][y*J>+i, ... ,x*J''''+^]^ 

> This corresponds to B{G — ij^) by (I30p 
Recompute i?~^ l> Using Sherman-Morrision- Woodbury formula 

Set k := k + fi > The number of copies of ij removed so far 

end if 

Set fi := fi/2 > The number of edges we try to remove is halved 

end while 
end for 
Return E > All removable copies of edges have been removed, so what remains is an /-factor 



This time we modify |17| Algorithm 1] in the same way as given by Algorithm [6j The time for 
updates remains the same, because a submatrix of size l^l x |5| is still updated in 0(|S'|'^) time. 
On the other hand, the cost we pay to find maximum number fj, of removable copies of edge e is 
0{{f{i) + /(j))"^)- The number of removed edges in the binary search forms a geometric series, so 
the cost is dominated by the first element, which in turn is smaller then the size of the submatrix 
to power of w. This time we obtain following bounds 

h{s) = Y,Ks-^) + 9{Si,S2) + o{\f{ST) 

i 
^^5^ /0(l/(^) + /(5)r) if |i?| = |5| = 1, 

' \Li,j9{R^,Sj) + 0{\f{R) + f{S)n otherwise. 

The solution for these equations gives an 0{(j)^) time bound. 

13 Finding perturbed factor weights 

This section shows how to compute the quantities w{Fy) and w{F'"), for all v (i V. (Recall from 
Section [7] that we are dealing with an /-critical graph; each v £ V has a maximum /^-factor F^ 
and a maximum /^'-factor F^'.) We start by considering simple graphs, and then comment on 
multigraphs. For simplicity assume that the weight function is non-negative, i.e., w : E ^ [0..W]. 
(If this is not the case, redefine w{ij) := w{ij) + W. This increases the weight of each /-factor by 
exactly Wf{V)/2. 
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Following (j25p . define 

'z'"i'^)xi^yi^ ij G E and (i, r, j, c) E P, 



B{G)^,r 



J,c 



otherwise, 



where z is a new indeterminate. For the next result we assume G has an /-factor. Theorem UT] 
shows that there is a mapping J^ from terms of det(i?(G)) onto $2(G). The degree of ^; in a term 
a equals the total weight of the edges used. This gives the following. 

Corollary 49. For a simple graph G that has an f -factor, deg^(det(i?(G))) is twice the weight of 
a maximum f -factor. 

Now suppose G is /-critical. For any v ^ V lei Gy be G with an additional vertex t joined 
to u G y by a zero weight edge. Set /(t) = 1. A maximum /-factor in G„ weighs the same as a 
maximum /^-factor in G. For the computation of /^-factors we need the following definition. Let 
G* be G with an additional vertex t that is connected to all vertices v £ V with zero weight edges. 
As previously, we set f{t) = 1. Let us denote by F* the maximum /-factor in G*. 

Lemma 50. deg2(adj(S(G*))t,,o,t,o) = w{F^) + w{Fy). 

Proof. Observe that 

adj(S(G,)),,o,t,o = (-l)"(*'0)+-(-,o) det(i?(G*)*'0'^'0), 

By Theorem 37] we know that det(-B(G*)) contains terms corresponding to elements of F2(G'^,). 
Hence, by the above equality, terms of adj(i?(G^,))t,^o,t,o correspond to elements -F2(G^,) that use 
edge tv, but with this edge removed. In other words terms of adj(i?(G*))„^o,t,o, are obtained by 
pairing an /-factor in G* and an /-factor in G^, and removing the edge tv. Similarly, as we 
observed in Corollary 1491 the degree of z encodes the total weight of elements of F2(G*). Moreover, 
the maximum elements are constructed by taking maximum /-factor of G* and maximum /-factor 
in Gy. As we already observed maximum /-factor in Gy is an maximum /^-factor in G^ so the 
theorem follows. ■ 

The above theorem leads to Algorithm [7] that computes weights of Fy, for all -y G y, in 0{Wn'^) 
time. 

Algorithm 7 An 0{W(j)^) time algorithm for finding weights of Fy, for all ?; G F, in a simple graph 

G+. 

Let B{Gif) he (f) X (p matrix representing G* 

Replace the variables in B{G^) for random elements from Zp for prime p = @{(j)'^) to obtain B 
Compute d := det{B) > requires 0(VF(/)'^) time using Theorem [3] 

{deg^{d) = deg2(det(i?(G*))) with probability > 1 — -ri) !> by Lemma [1] 

Set w{F^) := deg^(d)/2 > by Corollary |49] 

Compute a := adj(i?)ei^o = det(i?)i?~^ej^o '> requires 0{W(l)^) time using Theorem [3] 

for i; G y do 

(deg^(o„) = deg2(adj(B(G*))t,,o,t,o) with probability > 1 - ^) > by Lemma[T] 

Set w{Fy) := deg^iay) — w{F^) t> by Lemma [50l 

end for > by union bound all w{Fy) are correct with probability > 1 — ^ 
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For the computation of /'^-factors we need to proceed in slightly modified way. We construct 
G"" from G+ by: 

• adding new vertices tu and a zero weight edges ut^, for every vertex u G T^, 

• adding a new vertex t and zero weight edge tt^. 

Moreover, we define f'{v) = f{v) + 1, f'{ty) = 1 and f'{t) = 1. Finally, G* is obtained from G~^ 
in similar way, but with the difference that t is connected to all vertices ty. Again we observe, that 
the weight of the maximum /'-factor in G^ is equal to the weight of the maximum /^'-factor in G~^ . 
This allows us to prove the following. 

Lemma 51. deg,(adj(S(G*))i„,o,t,o) = w{F*)+w{F-). 

Proof. By Theorem [17] we know that det{B{G*)) contains terms corresponding to elements of 
F2{G*). Hence, terms of adj(i?(G*))j^^o,t,o correspond to elements F2{G*) that use edge tty, but 
with this edge removed. In other words terms of adi{B(G*))t^fl^v,o are obtained by pairing an /'- 
factor in G* and an /'-factor in G^ . As previously, the maximum weights elements are constructed 
by taking the maximum /'-factor of G* and the maximum /'-factor in G^ . As observed above 
/'-factors correspond to /''-factors in G"*". ■ 

This leads to 0{Wn'^) time algorithm for computing weights of F"" , for all v £ V. 

Algorithm 8 An 0{W(j)''^) time algorithm for finding weights of F"" , for all v G V, for a simple 
graph G"^. 

1 

2 

3 



5: 
6: 

7: 

8: 

9: 
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Let B{G*) he (j) X (p matrix representing G* 

Replace the variables in B{G*) for random elements from Zp for prime p = 0((/)^) to obtain B 
Compute d := det(i?) I> requires 0{W(f)^) time using Theorem [3] 

(degj,(d) = deg2(det(i?(G*))) with probability > 1 - ^j) i> by Lemma [1] 

Set w{F*) := deg^(d)/2 > by Corollary |49] 

Compute a := adj(i?)ei^o = det(i?)i?~^ej^o l> requires 0{W(j)^) time using Theorem [3] 

for f G y do 

(deg^(afj = deg^(adj(5(G*))t„,o,i,o) with probability > 1 - jj) > by Lemma [U 

Set w^F"") := deg;,{ay) - w{F*) > by Lemma EH 

end for > by union bound all w{F^') are correct with probability > 1 — ^ 



13.1 Multigraphs 

Let G = {V, E) be a multigraph, and let w : i? x /c — )■ Z be the edge weight function. This function 
assigns weight w;(e, k) to the /c'th copy of e G i?. Joining ideas from (j34p and (f26|) we define 

{E0? z'"^^i^^)3^J'^yf^ ij G E and (i, r, j, c) G P, 
- E0? z-'^^^^a^^'^vV^'' ij G E and (j, c, i, r) G P, 
otherwise, 

where z is a new indeterminate. Theorem |47] shows that there is a mapping F from terms of 
det (i?(G)) onto F2{G). Observe that the construction from previous section requires only the 
existence of such mapping, so it can be used for multigraphs as well. 
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14 Conclusions and open problems 

This paper presents new algebraic algorithms for the fundamental problems of 6-matching, undi- 
rected single-source shortest paths, and /-factors. Some intriguing open problem and challenges 
emerge from this study: 

• The matrices we construct for unweighted /-factors have a very special block structure. Can 
this block structure be exploited to obtain faster algorithms, e.g., time 0{(j)^~^7i)7 

• Can the running time of our algebraic max-flow algorithm be improved, perhaps by combining 
it with scaling techniques? Can scaling be used in the non-bipartite algorithms? 

• We gave first algebraic algorithms for simple 2- factors. Are there algebraic formulations for 
triangle-free or square-free 2-factors? If so one expects the resulting algorithms to be simpler than 
existing combinatoric ones. 

• What is the complexity of all-pairs undirected shortest distances on conservative graphs? Can 
0{Wn'^) time be achieved, as in the case of non-negative weights? 
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A Allowed edges 



b) 



Figure 8: The only 2- factor is shown on panel a), whereas the only /^^^-factor is shown on panel b). 
Edge uv is not allowed in any 2-factor, although /^^^-factor does exist. 

Observe that for any edge ij as long as f{i) = f{j) = 1 the existence of /jj-factor is equivalent to 
the fact that ij is allowed. However, as shown on Figure [U the edge uv is not allowed, although the 
fu^v-^actor does exit. The reason for this is that we are trying to use edge uv twice. Hence, similar 
criteria as in Corollary [9] does not hold in non-bipartite case. 
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